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The  research  program  described  in  this  ELnal  Report  has  included  work  on  boundary 
layer  studies,  current  distribution  and  current  leakage  in  MHD  generators,  and  boundary 
layer  diagnostics.  A  computer  program  which  has  been  used  successfully  to  calculate 
conventional  turbulent  boundary  layers  has  been  modified  by  adding  the  MHD  terms  in  tht 
momentum  and  energy  equations .  ^ This  modified  program  has  been  coupled  to  a  one-dimen¬ 
sional  core-flow  program  to  give  boundary  layer  parameters  in  an  MHD  generator.  Exper¬ 
iments  were  made  in  arc-heated  argon  with  small  additions  of  N2  and  CO2  to  delineate 
the  effect  of  these  components  on  electrode  voltage  drops.  The  effect  was  larger  than 
could  be  predicted  from  ambipolar  diffusion  theory.  Analysis  and  experiments  made  wit! 
electrodes  of  several  sizes  showed  the  important  influence  of  boundary  layer  and  elec¬ 
trode  temperatures.  Because  the  boundary  layer  temperature  is  coupled  to  the  electrodt 
temperature,  small  cold  electrodes  performed  better  than  large  cold  electrodes  in  the 
Stanford  MHD  generator. 

.Current  patterns  in  conducting  wall  generators  were  measured  and  showed  good  agree 
ment  with  theoretical  predictions  for  the  anode.  The  surface-sheath  voltage  drop  was 
responsible  for  less  good  agreement  at  the  '’athode.  It  was  shown  that  a  large  fractior 
of  the  current  entered  the  sidewalls  of  a  conducting  wall  generator.  Experiments  also 
indicated  that  Faraday  electrodes  which  included  part  of  the  sidewall  gave  a  greater 
power  output  than  conventional  flat  electrodes.  Axial  current  leakage  in  the  Hall 
direction  was  calculated  and  measured  in  a  simulated  MHD  generator  using  argon  with  up 
to  20%  N2  diluent.  Agreement  between  experiment  and  theory  was  satisfactory.  For  com¬ 
bustion  gases  with  diffuse  currents  the  theory  predicts  larger  Hall  voltages  before 
breakdown  than  have  been  observed  in  practice.  The  Hall  voltage  limitation  is  believe 
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associated  with  arcing  phenomena. 

I  A  new  method  which  uses  Faraday  splitting  of 
laser  modes  is  suggested  for  the  measurement  of 
electron  concentration  in  the  boundary  layer.  A 
laser  doppler  velocimeter  and  its  applications  to 
boundary  layer  velocity  measurements  in  the 
Stanford  generator  is  described. 
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FOREWORD 


The  research  described  in  the  Final  Report  was  performed  in  the  High 
Temperature  Gasdynamics  Laboratory  of  Stanford  University,  Stanford,  California, 
under  contract  AF  F33615-69-C-1171,  project  number  3145,  task  number  314526. 

The  program  was  supported  by  the  Air  Force  Aero  Propulsion  Laboratory,  Wright 
Patterson  Air  Force  Base,  Ohio.  Mr.  Robert  F.  Cooper  and  Mr.  Paul  Lindquist, 
POP-2,  served  as  technical  managers.  This  report  covers  the  period  9/15/68 
to  7/31/71  and  was  submitted  on  January  24,  1972. 

The  authors  would  like  to  point  out  the  very  important  contributions  made 
by  Stanford  graduate  students  during  this  research  program.  In  particular  Mr. 

S.  J.  Morris  was  responsible  for  the  computer  program  described  in  Section  2.1, 
Mr.  G.  A.  DuBro  made  the  measurements  of  Section  2.2,  and  Dr.  E.  S.  Rubin  per¬ 
formed  the  calculations  and  experiments  presented  in  Section  2.3.  Mr.  R.  M 
Cima  and  Mr.  K.  E.  Berry  assisted  with  the  current  distribution  work  described 
in  Section  3.1.  Dr.  J.  C.  Cutting  was  responsible  for  the  current  leakage 
program  discussed  in  Section  3.2.  The  boundary  layer  diagnostic  work  given  in 
Sections  4.1  and  4.2  was  conducted  by  Dr.  S.  A.  Self,  Senior  Research  Associate. 
Dr.  Robert  Kessler  assisted  with  much  of  the  laboratory  program  and  Mr.  Frank 
Levy's  assistance  in  the  experimental  work  was  invaluable. 


ABSTRACT 


The  research  program  described  in  this  Final  Report  has  included  work 
on  boundary  layer  studies,  current  distribution,  and  current  leakage  in  magneto¬ 
hydrodynamic  (MHD)  generators,  and  boundary  layer  diagnostics.  A  computer 
program  which  has  been  used  successfully  to  calculate  conventional  turbulent 
boundary  layers  has  been  modified  by  adding  th~  MHD  terms  in  the  momentum  and 
energy  equations.  This  modified  program  has  been  .oupled  to  a  one-dimensional 
core-flow  program  to  give  boundary  layer  parameters  in  an  MHD  generator.  Exper¬ 
iments  were  made  in  arc-heated  argon  with  small  additions  of  and  to  cselin 
eate  the  effect  of  these  components  on  electrode  voltage  drops.  The  effect 
was  larger  than  could  be  predicted  from  ambipolar  diffusion  theory.  Analysis 
and  experiments  made  with  electrodes  of  several  sizes  showed  the  important 
influence  of  boundary  layer  and  electrode  temperatures.  Because  the  boundary 
layer  temperature  is  coupled  to  the  electrode  temperature,  small  cold  elec¬ 
trodes  performed  better  than  large  cold  electrodes  in  the  Stanford  MHD  generator 
Current  patterns  in  conducting  wall  generators  were  measured  and  showed 
good  agreement  with  theoretical  predictions  for  the  anode.  The  si  .face-sheath 
voltage  drop  was  responsible  for  less  good  agreement  at  the  cathode.  It  was 
shown  that  a  large  fraction  of  the  current  entered  the  sidewalls  of  a  conducting 
wall  generator.  Experiments  also  indicated  that  Faraday  elect\odes  which 
included  part  of  the  sidewall  gave  a  greater  power  output  than  conventional  flat 
electrodes.  Axial  current  leakage  in  the  hall  direction  was  calculated  and 
measured  in  a  simulated  MHD  generator  using  argon  with  up  to  20%  ^  diluent. 
Agreement  between  experiment  and  theory  was  satisfactory.  For  combustion  gases 
with  diffuse  currents  the  theory  predicts  larger  Hall  voltages  before  breakdown 
than  have  been  observed  in  parctice.  The  Hall  voltage  limitation  is  believed 
associated  with  arcing  phenomena. 
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A  new  method  which  uses  Faraday  splitting  of  laser  modes  is  suggested 
for  the  measurement  of  electron  concentration  ii\  the  boundary  layer.  A  laser 
doppler  velocimeter  and  its  applications  to  boundary  layer  velocity  measurements 
in  the  Stanford  generator  is  described. 
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1.0  INTRODUCTION 


The  High  Temperature  Gasdynamics  Laboratory  of  Stanford  University  has 
been  investigating  losses  originating  near  or  in  the  boundary  layer  of  an 
Magnetohydrodynamic  (MHD)  generator  with  a  view  to  reducing  these  losses.  The 
research  program  has  been  directed  at  experimental  and  analytical  studies  of 
the  electrode  boundary  layer,  of  the  current  distribution  and  leakage  as  influ¬ 
enced  by  the  boundary  layers,  and  at  diagnostic  techniques  which  appear  capable 
of  the  necessary  spatial  resolution  to  be  useful  in  the  boundary  layer. 

Certain  parts  of  the  effort  have  been  completed  and  have  been  reported 
ir>  Air  F: rce  Aero  Propulsion  Laboratory  reports .  Other  parts  of  smaller  scope 
have  bf  an  reported  in  the  technical  literature.  Work  which  was  underway 
near  the  conclusion  of  the  contract  has  not  been  previously  reported  and  is 
given  in  more  detail  in  this  final  report. 

The  following  publications  have  resulted  from  the  program  conducted  under 
contract  F33615-72-C-1171: 

"Effects  of  Electrode  Size  on  the  Performance  of  a  Combustion-Driven  MHD 
Generator,"  E.  S.  Rubin,  SU-IPR  Report  #331,,  August  1969. 

"A  Computer  Program  for  the  Calculation  of  the  Electron  Transport  Coefficients," 
M.  Mitchner,  SU-IPR  Report  #280,  August  1969. 

"Axial  Current  Leakage  in  Segmented  MHD  Generators,"  J.  C.  Cutting,  SU-IPR 
Report  #406,  January  1971. 

"Measurement  of  Current  Distributions  and  the  Effect  of  Electrode  Config¬ 
uration  on  MHD  Generator  Performance,"  R.  H.  Eustis  and  R.  Kesslet,  Fifth 
Inti.  Conf.  on  MHD  Electrical  Power  Generation,  I_,  April  1971. 

"Axial  Current  Leakage  in  Segmented  MHD  Generators,"  J.  C.  Cutting  and  R.  H. 
Eustis,  Fifth  Inti.  Conf.  on  MHD  Electrical  Power  Generation,  _1,  April  1971. 

"Effects  of  Electrode  Size  on  the  Performance,  of  a  Combustion-Driven  MHD 
Generator,"  E.  S.  Rubin  and  R.  H.  Eustis,  AIAA  Journal,  Vol.  9,  No.  6, 

June  1971. 

"Axial  Current  Leakage  Between  Electrodes  in  MHD  Generators,"  J.  C.  Cutting, 

R.  H.  Eustis,  Eleventh  Symposium  on  Engineering  Aspects  of  Magnetohvdrody- 
namics,  1970. 
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"Electrode  Size  Effects  in  Combustion-Driven  MHD  Generators,"  E.  S.  Rubin, 

R.  H.  Eustis,  Eleventh  Symposium  on  Engineering  Aspects  of  Magnetohydrody¬ 
namics,  1970 

"Effects  of  Electrode  Size  on  MHD  Generator  Performance,"  E.  S,  Rubin,  R.  H. 
Eustis  Tenth  Symposium  on  Engineering  Aspects  of  Magnetogydrodynamics ,  1969. 

"Current  Distribution  in  Conducting  Wall  MHD  Generators,"  R.  H.  Eustis,  R.  M. 
Cima,  and  K.  E.  Berry,  Eleventh  Symposium  on  Engineering  Aspects  of  Magne¬ 
tohydrodynamics,  1970. 


The  present  report  which  constitutes  the  final  report  for  contract  F33615- 
69-C-1171  is  divided  into  three  major  sections.  The  first  deals  with  the  boundary 
layer  analysis  and  with  experiments  related  to  effects  influencing  the  near¬ 
electrode  region.  The  second  section  is  concerned  with  the  calculation  and 
measurement  of  current  distribution  in  a  conducting  wall  generator  and  with 
current  leakage  in  the  Hall  direction.  The  last  section  deals  with  diagnostic 
techniques.  A  brief  summary  concludes  the  report. 


2.0  BOUNDARY  LAYER  STUDIES 

In  this  section  a  description  of  the  analysis  of  the  turbulent  MHD 
boundary  layer  and  some  computer  results  are  given,  the  measurements  relating 
to  the  effect  on  electrode  voltage  drops  due  to  molecular  species  are  described, 
and  the  results  of  tests  with  different  electrode  sizes  are  presented. 

2.1  Computer  calculations  of  boundary  layer  parameters 

In  order  to  take  advantage  of  the  large  amount  of  work  which  has  been 
done  in  developing  computer  programs  to  calculate  conventional  turbulent  boundary 
layers,  the  excellent  program  of  Patankar  and  Spalding  [l]*as  modified  exten¬ 
sively  by  Kays  and  coworkers  [2  ]  has  been  used  as  a  basis  for  our  calculations. 
The  basic  equations  have  been  modified  to  include  the  MHD  terms  as  follows: 


momentum  equation 

3u  3u, 

(ut—  +  vr— )  -  j 
dx  dy  y 


M  +  2~ 

9x  9x 


,  3u. 

V 


y  -  momentum  equation 


*Numbers  in  square  brackets  refer  to  the  References. 
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energy  equation 
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3h 
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pu 


where 
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+  pv 
=V 
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3u 

xy  -h*eff  3y 
3T 


3  -»•-*■ 

t—  (uT  )  +  E  •  j  +  tt- 
3y  xy  J  3y 


eff  3y 


For  computational  convenience  the  above  equations  have  been  rewritten  in  terms 
of  the  variables 


ip  m  -pv  and  ip  =  pu 
X  y 


as  suggested  by  von  Mises  [1]  in  order  to  automatically  satisfy  the  continuity 
equation.  The  Kays  -  Patankar  -  Spalding  computer  program  has  been  altered  to 
account  for  the  extra  terms  due  to  the  MHD  effects  and  runs  have  been  made  as 
described  later. 

It  is  necessary  to  give  the  computer  program  a  pressure  distribution 
which  is  derived  from  a  one-dimensional  analysis  of  the  core  flow.  This 
analysis  uses  the  boundary  layer  displacement  thickness  as  a  correction  for 
the  core  flow  area  distribution  so  iteraction  between  the  core  flow  and  bound¬ 
ary  layer  programs  is  required.  The  core  flow  analysis  makes  use  of  the 


standard  momentum,  energy  and  .. ontinuity  equations 

pu  ~  ^  -  (J  x  t) 

K  dx  dx  VJ  'x 


dh 


pu 


dx 


=  E  •  j  +  q'" 

J  Mrad 


>Aeff  V 


constant 


wher~  Is  the  flow  area  corrected  for  the  boundary  layer  and  the  radiation 


term  includes  band  radiation  from  and  HjO  and  line  radiation  from  the 


potassium  seed. 

An  erperirac^tal  program  has  been  planned  to  measure  the  boundary  layer 
velocity  and  temperature  profiles  in  an  MHD  generator  and  the  proposed  oper¬ 
ating  conditions  are  shown  in  Figure  1.  An  illustration  of  the  application 
of  the  computer  program,  the  calculated  profiles  and  other  parameters  for 
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the  proposed  experiment  are  shown  in  Figures  2  to  4.  It  is  seen  that  two  major 
test  conditions  are  planned  —  one  subsonic  and  the  other  supersonic.  The 
apparatus  consists  of  the  burner  and  plenum  followed  by  a  run-in  section  0.6 
meter  long  and  the  test  section  0.7  meter  long  which  is  inserted  in  the  magnet. 
The  run-in  section  is  to  thicken  the  boundary  layer  to  facilitate  measurements 
and  to  more  closely  simulate  larger  scale  generators. 

The  variation  of  the  core  velocity  and  the  pressure  with  axial  distance 
measured  from  the  nozzle  is  shown  in  Figure  2.  The  dotted  lines  which  show 
the  case  for  B  *  2  tesla  indicate  opposite  effects  for  subsonic  and  supersonic 
flow  as  is  expected.  Velocity  profiles  near  the  end  of  the  test  section  are 
shown  in  Figure  3  and  are  plotted  in  a  normalized  fashion  in  Figure  4.  It 
is  clear  that  the  profile  is  modified  ..  shape  as  well  as  in  the  magnitude 
of  velocities  by  the  MUD  effect. 

Measurements  of  the  velocity  near  the  outer  edge  of  the  boundary  layer 
were  conducted  under  contract  F33635-69-C-1171  as  discussed  in  section  4.2  of 
this  report.  Further  measurements  are  being  conducted  under  contract  F33615- 
72-C1088. 

2.2  Investigation  of  the  effpcts  of  molecular  gases  on  electrode  voltage 
losses . 

A  series  of  experiments  designed  to  obtain  an  understanding  of  the  nature 
of  electrode  voltage  drops  in  combustion-gas  MBD  generation  were  performed.  'Hie 
objective  of  these  experiments  and  associated  theoretical  work  was  to  be  able  to 
predict  the  magnitude  of  this  important  loss  mechanism  for  proposed  MHD  generator 
designs. 

A  series  of  three  experiments  have  shown  that  the  addition  of  small  amounts 
of  carbon  dioxide  into  an  argon-potassium  seeded  plasma  has  substantial  effect 
on  electrode  voltage  drop  behavior.  With  the  addition  of  5%  CO^,  cathode  voltage 
drops  comparable  to  those  exhibited  in  combustion-gas  experiments,  have  been 
measured.  (See  Figure  5  ) 
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Figure  5.  Cathode  Voltage  Drop  with  Argon  Plus  Additives 
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The  experimental  apparatus  consisted  of  an  argon  arc  jet,  a  plenum  chamber 
where  controlled  amounts  of  potassium  seed  and  molecular  gases  were  introduced 
and  a  quartz  test  section  in  which  a  discharge  was  imposed  across  multiple 
pairs  of  segmented  electrodes.  Voltage  drops  were  measured  by  means  of  an 
array  of  probes  located  between  the  discharge  electrodes.  The  gas  temperature 
was  in  the  range  1900°K  and  the  electrode  temperature  was  approximately  1400°K. 

Potassium  seeding  was  0.2Z  (mass  fraction)  and  CO2  addition  was  varied  from  OX 

2 

to  51.  With  externally  applied  voltages,  current  densities  up  to  0.5  amps/cm 
were  measured. 

Figure  5  shows  the  relative  magnitudes  of  electrode  voltage  drops  in  various 
experimental  situations.  The  combustion  gas  data  is  taken  from  Kessler  [3  ]  and 
the  nitrogen-argon  data  was  provided  by  Cutting  [4].  Figures  6and7show  the 
results  of  cathode  and  anode  voltage  drops  respectively  versus  current  density 
for  the  first  two  experiments.  Similar  behavior  was  obtained  in  the  third 
experiment,  but  different  gas  conditions  did  not  allow  comparison  on  these  figures. 

From  these  data  the  following  observations  may  be  made: 

(1)  For  specified  gas  conditions,  i.e.,  for  specified  electrode  and  gas 
temperatures  and  seed  fraction,  increasing  addition  of  CO^  resulted 
in  higher  voltage  drops . 

(2)  Electrode  voltage  drops  were  linear  with  current  density  for  small 
current  densities. 

(3)  Relatively  small  additions  of  CO^  have  substantial  effects  on  electrode 
behavior,  especially  at  the  cathode,  (e.g.,  .01%  CO^  addition  essen¬ 
tially  doubles  the  cathode  voltage  drop.) 

(4)  Cathode  voltage  drops  are  several  times  larger  than  the  correspond  mg 
anode  drops. 

A  theoretical  interpretation  of  these  results  has  been  attempted  on  the 
basis  cr  an  ambipolar  -  diffusion  model  which  accounts  for  the  boundary-layer 
temperature  profile.  This  model  has  been  moderately  successful  in  predicting 
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Figure  6.  Effect  of  CO  on  Cathode  Voltage  Drop 
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ELECTRODE  VOLTAGE  DROP  ON  ANODE 


the  observed  behavior,  but  it  appears  that  additional  features  Bust  be  incorpor¬ 
ated  to  obtain  an  accurate  description. 

This  work  will  be  completely  describe*,  in  a  forthcoming  Stanford  report 
which  will  be  sent  to  the  Aero  Propulsion  Laboratory. 

2.3  Measurement  of  the  effect  of  electrode  size  on  the  performance  of  an 
MUD  generator 

Experiments  were  conducted  to  investigate  the  effects  of  electrode  size  on 
the  performance  of  a  combustion-driven  MHD  generator.  Electrode  sizes  with 
length-to-pitch  ratios  of  0.23,  0.50,  and  0.79  were  tested  in  a  generator  section 
simulated  by  three  molybdenum  electrode  pairs  located  at  the  downstream  end  of 
an  MHD  channel.  Voltage  probe  data  were  obtained  for  electrode  surface  temper¬ 
atures  between  500°K  and  1600“K.  with  gas  conductivity  in  the  electrode  boundary 
layer  regions  established  independently  by  either  a  hot  (2200°K  MgO  brick)  or 
cold  (750°K  water-cooled  plate)  upstream  wall.  At  a  given  load  current,  larger 
electrodes  were  found  to  have  lower  voltage  losses  for  similar  conditions  of 
surface  and  gas  boundary- layer  temperature.  For  dissimilar  boundary  layers, 
reflecting  coupling  between  electrode  and  boundary-layer  temperatures  in  a  cold- 
electrode  generator,  total  voltage  losses  for  a  large  electrode  pair  were  equal 
to  or  greater  than  those  of  a  small  electrode  pair  at  the  same  surface  tempera¬ 
ture.  Experimental  boundary- layer  resistances  at  anodes  were  in  good  agreement 
with  analytical  predictions  based  on  both  an  approximate  and  a  more  exact  gas- 
dynamic  model,  for  electrode  temperatures  above  about  800°K.  At  lower  temperatures 
anodes  were  believed  to  be  influenced  by  solidification  of  seed  compounds  on 
the  electrode  surfaces,  resulting  in  higher  losses.  For  cathodes,  the  boundary- 
layer  resistance  was  gasdynamic  for  the  case  of  a  hot  upstream  wall,  for  electrode 
current  densities  below  that  of  saturated  thermionic  emission  at  the  electrode 
temperature.  In  all  other  cases,  the  presence  of  cathode  surface-sheath  effects 
was  indicated. 
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This  work  has  been  reported  in  detail  in  Air  Force  Technical  Report  No. 
AFAPL-TR-69-72  and  in  AIAA  Journal,  Vol.  9,  No.  6,  June  1971,  pp.  1162-29. 

3.Q  CURRENT  DISTRIBUTION  AND  LEAKAGE  STUDIES 

Two  major  efforts  of  contract  F33615-69-C-1171  were  devoted  to  the  calcula¬ 
tion  and  measurement  of  electrode  current  distribution  in  a  conducting  wall 
generator  and  to  the  calculation  and  measurement  of  Hall  current  leakage  in  a 
simulated  MHD  generator.  Both  of  these  programs  have  been  reported  in  detail 
in  the  literature  and  are  summarized  below. 

3.1  Electrode  current  distribution 

An  experimental  program  was  conducted  to  investigate  the  jerforaance  of 
a  combustion-driven  MHD  generator  with  conducting  sidewalls  and  with  various 
electrode  configurations.  The  current  distribution  to  the  electrodes  and  conducting 
sidewalls  of  a  Hall  generator  was  measured.  The  Hall  generator  was  operated 
with  both  conducting  and  with  insulating  sidewalls;  no  significant  performance 
differences  were  found.  Significant  current  was  found  to  flow  to  the  sidewall 
segments  of  the  conducting-wall  Hall  generator,  especially  to  those  segments 
immediately  adjacent  to  the  normal  electrodes.  An  experiment  was  made  with  a 
Faraday  generator  in  which  the  electrodes  included  the  immediately  adjacent 
sidewall  segments.  This  generator  produced  about  25%  more  power  than  a  conven¬ 
tional  Faraday  generator  with  flat  electrodes. 

This  type  of  electrode  in  which  part  of  the  current  enters  (or  leaves)  the 
electrode  parallel  to  the  magnetic  field  may  be  of  special  interest  for  high 
magnetic  field  generators.  Current  concentrations  should  be  reduced  over  those 
on  flat  electrodes  end  higher  Hall  voltages  may  therefore  be  possible. 

A  computation  procedure  was  developed  for  calculating  the  current  distri¬ 
bution  in  the  plane  perpendicular  to  the  velocity  vector.  Appropriate  turbulent 
velocity  and  temperature  profiles  were  assumed  and  the  results  agreed  well  with 
the  measured  anode  current  distribution.  The  agreement  of  measured  and  predicted 
current  distributions  at  the  cathode  was  less  satisfactory  probably  because  of 
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the  cathode  surface-sheath  voltage  drops  which  were  not  accounted  for  in  the 
model. 

This  work  has  been  reported  in  detail  in  the  following  publications: 

Eustis,  Cima,  Berry,  "Current  Distribution  in  Conducting  Wall  Generators," 
Eleventh  Symposium  on  Engineering  Aspects  of  Magnetohydrodynamics,  1970,  and 
Eustis,  Kessler,  "Measurement  of  Current  Distributions  and  the  Effect  of  Electrode 
Configuration  on  MHD  Generator  Performance,"  Fifth  Inti.  Conf.  on  MHD  Elec¬ 
trical  Power  Generation,  Xi April  1971. 

3.2  Axial  current  leakage  in  segmented  MHD  generators 

Experiments  were  performed  in  a  simulated  linear,  segmented  MHD  generator 
to  investigate  the  phenomenon  of  axial  current  leakage  between  electrodes.  An 
argon-potassium  plasma  was  used  with  the  addition  of  up  to  20  percent  nitrogen 
to  suppress  non-equilibrium  effects.  Electric  and  magnetic  fields  were  applied 
to  provide  a  maximum  average  transverse  current  density  of  2  amperes  per  sq  cm 
and  a  maximum  Hall  parameter  of  6. 

The  generator  was  composed  of  six  tungsten  electrode  pairs  inserted  into 
a  rectangular  quartz  channel.  The  measuring  or  "test"  electrode  was  electrically 
guarded  and  relaxation  effects  minimized  by  three  upstream  electrode  pairs  and 
two  pairs  downstream.  The  anode  and  cathode  of  this  "test"  electrode  were 
divided  into  three  segments  each  and  the  current  direction  and  magnitude  measured 
at  eacn  segment. 

Axial  current  leakage  between  adjacent  electrodes  was  measured  and  found 
to  be  a  function  of  increasing  transverse  current  density.  The  "critical" 
transverse  current  density  necessary  to  cause  leakage  decreased  with  increased 
Hall  parameter,  decreased  plasma  temperature  and  decreased  percentage  of  nitro¬ 
gen  additive.  Leakage  currents  in  excess  of  50%  of  transverse  currents  were 
measured  with  no  indication  of  breakdown  or  arcing  between  adjacent  electrodes. 
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Hall  voltages  were  observed  to  have  an  asymptotic  limit  with  increasing  transverse 
current  density.  The  value  of  this  limit  was  a  function  of  nitrogen  percentage 
and  Hall  parameter. 

A  computer  program  for  calculation  of  electrical  current  distributions 
in  a  segmented  MHD  device  was  written  to  interpret  and  predict  the  experimental 
data.  The  program  consisted  of  a  coupled  solution  of  the  current  stream  function 
equations  and  the  electron  energy  and  continuity  equations.  It  was  found  that 
the  Hall  voltage  data  and  the  measured  cathode  current  phenomena  were  well 
predicted  by  the  theory.  Prediction  of  current  leakage  at  the  anode  wall  did 
not  compare  as  well  with  the  measured  results.  The  calculations  indicated  that 
the  increase  in  electrical  conductivity  along  the  electrode  walls  was  responsible 
for  axial  current  leakage. 

For  combustion  gases  where  the  electrical  conductivity  is  nearly  independent 
of  current  density  (for  modest  densities)  the  calculations  showed  that  current 
leakage  would  not  occur  for  Hall  fields  several  times  those  achieved  in  practice. 
It  is  believed  that  the  Hall  breakdown  which  has  occurred  in  generators  is  asso¬ 
ciated  with  arcing  rather  than  a  limitation  which  would  occur  for  diffuse 
currents. 

This  work  has  been  reported  completely  in  Air  Force  Technical  Report  No. 
AFAPL-TR-71-23,  January,  1971,  by  John  C.  Cutting  and  in  the  paper  by  Cutting 
and  Eustis,  "Axial  Current  Leakage  Between  Electrodes  in  MHD  Generators," 

Eleventh  Symposium  on  Engineering  Aspects  of  Magnetohydrodynamics,  1970. 

4.0  MHD  GENERATOR  BOUNDARY  LAYER  DIAGNOSTICS 

In  an  MHD  generator  boundary  layer  three  of  the  important  quantiti 
which  need  experimental  determination  are  gas  temperature,  gas  velocity,  and 
electron  number  density.  The  measurement  of  gas  temperature  was  described  in 
a  report  for  contract  AF  F33615-67-C-1127 ,  Techncial  Report  No.  AFAPL-TR-68-141 . 

In  the  sections  below  a  suggested  new  method  of  measuring  electron  number  density 
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with  good  spatial  resolution  is  descr  e  and  the  development  and  testing  of 

a  laser  doppler  velocimeter  is  given.  This  latter  topic  was  a  major  diagnostic 

effort  of  contract  AF  F33615-69-C-1171. 

4.1  Electron  number  density  measurement 

A  new  method  is  described  for  measuring  electron  concentration  in  a 
magnetoplasma  by  observing  the  beat  frequency  due  to  Faraday  splitting  of 
a  laser  mode  when  the  plasma  is  situated  in  the  laser  cavity  with  the  magnetic 
field  along  the  laser  axis.  It  appears  to  offer  a  number  of  advantages  over  con¬ 
ventional  laser  interferometric  methods  with  regard  to  sensitivity  and  spatial 
resolution. 

In  conventional  interferometric  methods  for  determining  electron  number 
density  n  at  microwave  or  optical  frequencies  [5],  the  minimum  density  (or  den¬ 
sity  change)  that  can  be  measured  is  set  by  the  requirement  that  the  density 

a 

line  integral  ndx  =  uZ  should  gi\e  a  detectable  phase  shift.  For  an  unmagne- 

2  2  1/2  1 /? 

tized  plasma  the  refractive  index  is  y  =  (1  -  f  /f  )  =  (1  -  n/n  )  where 

P  c 
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fp  =  8.97  x  10  n  is  the  plasma  frequency  and  n£  *  1.24  x  10  f  is  the  cut¬ 
off  density  for  the  source  frequency  x  =  c/X  (cgs  units).  For  f  »  f^  the  phase 
shift  is  given  by 
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and  the  minimum  detectable  density  line  integral,  or  resolution  is 
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where  (A<J>/2Tr)ro^n  is  the  minimum  detectable  fractional  fringe  shift.  In  practice 
this  is  about  (1/10)  fringe,  though  it  may  b*  improved  to  ~ (1/ 100)  fringe  by  various 
sophisticated  techniques. 

For  laboratory  plasmas  of  dimensions  typically  Z  =  1  -  100  cm,  it  transpires 

that  microwaves  (up  to  the  maximum  practical  frequency  -10^  Hz)  are  satisfactory 

13  -3 

for  densities  n  <  10  cm  ,  while  visible  or  near  infrared  lasers  are  useful  for 


densities  n  >  10^  cm  ^ . 


For  the  important  intermediate  range. of  densities 
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10  <  n  <  10  ,  of  particular  interest  for  fusion  and  MHD  studies,  far  infra¬ 

red  lasers  (e.g.  HCN  at  337  ym)  have  recently  been  employed  [6-8  ]. 

While  resolution  in  density  is  the  primary  criterion  in  the  choice  of  fre¬ 
quency,  spatial  resolution  is  also  important.  The  narrow,  low-divergence  beams 
of  visible  and  near  infrared  lasers  give  good  spatial  resolution  but  are  insuf¬ 
ficiently  sensitive  except  for  very  high  densities.  It  is  clear  that  if  the 
sensitivity  could  be  markedly  increased,  then  visible  or  near  infrared  lasers 
could  be  used  in  the  intermediate  range  of  densities,  with  a  great  gain  in 
spatial  resolution.  Moreover,  visible  and  near  infrared  lasers,  together  with 
the  associated  optical  techniques,  are  more  highly  developed  and  readily 
available  than  their  far  infrared  counterparts. 

A  very  great  increase  in  sensitivity  may  be  achieved  [5,9]  by  incorpora¬ 
ting  the  plasma  into  the  laser  cavity  so  that  it  controls  the  frequency.  For 
a  mirror  spacing  L,  the  frequency  shift  Af  produced  by  the  plasma  is  given  by 

(3) 

while  the  resolution  is 

~  3.7  *  102  (f)  4f„ln  (4) 

where  q  ■  (2L/A  is  the  axial  mode  number.  The  minimum  frequency  shift  Af^^ 
that  can  be  detected  is  determined  by  the  frequency  stability  of  the  laser  and 
of  the  frequency  reference,  which  may  be  either  a  cavity  resonator  (Fabiy-Perot 
etalon)  or  another  identical  laser.  In  practice  far  greater  sensitivity  is 
obtained  by  heterodyning  against  a  second  laser. 

By  using  two  He-Ne  lasers,  Johnson  [5,9]  has  been  able  to  measure  fre¬ 
quency  shifts  on  repetitive  signals  (pulsed  plasma)  as  low  as  1  kHz,  corres¬ 
ponding  to  a  density  of  2  x  10^  cm  ^  in  a  plasma  30  cm  long.  To  achieve  this 
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resolution  interfarometrically  would  require  detection  of  a  fringe  shift  of 
10  6.  However,  this  great  sensitivity  is  achieved  only  at  the  expense  of 

f 

severe  difficulties  associated  with  thermal  and  mechanical  instability  of 
the  lasers.  Furthermore,  because  under  these  conditions  the  free  electron 
contribution  to  the  refractive  index  is  less  than  IQ-1®,  it  may  be  completely 
masked  by  the  contribution  due  to  excited  atoms  and  ions.  Johnson  overcame 
this  problem  by  simultaneous  measurements  at  A  ±  0.633  and  1.152  ym,  so  that 
by  subtraction  the  frequency-dependent  part  of  the  refractive  index  due  to 
free  electrons  was  obtained. 

I 

In  the  new  method  proposed  here,  we  take  advantage  of  the  increased 

:  i  ■ 

sensitivity  resulting  from  placing  the  plasma  in  the  laser  cavity  but,  by 
using  Faraday  aplitting,  we  avoid  the  complication- of  a  second  laser  and  the 
associated  stability  problems. 

For  propagation  along  the  direction  of  a  static  magnetic  field  B,  the 
plasma  refractive  index  is 


where  f  is  the  electron  cyclotron  frequency,  and  the  signs  correspond  to 
the  left  (+)  and  right  (-)  circularly  polarized  waves  defined  relative  to  the 

i 

direction  of  B. 

In  the  absence  of  a  gyrotropic  medium  a  given  laser  mode  TEM^^  is  linearly 
polarized  and  may  be  regarded  as  the  superposition  of  right  and  left  circularly 
polarizea  modes.  The  inclusion  of  an  axially  magnetized  plasma  in  the  laser 
cavity  is  just  Af  given  by  Eq.  3,  while  the  differential  shift  if  is  given  by 
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Thus  the  differential  shift  is  smaller  than  the  average  shift  by  a  factor  (2fc/f), 
while  the  resolution  is 

«  3-T  X  1q2(i)(t)  •  (7) 

Combining  the  two  component  modes  it  is  seen  that  they  correspond  to  a  composite 
plane  polarized  mode  whose  plane  of  polarization  rotates  at  a  frequency  (6f/2). 
Hence,  by  passing  the  laser  output  through  a  linear  polarizer,  the  frequency 
difference  may  be  detected  directly  as  an  amplitude  modulation  at  frequency  6f. 

Thermal  and  mechanical  fluctuations  will  affect  the  component  modes  equally, 
so  that  the  problems  of  frequency  stability  associated  with  the  beating  of  two 
lasers  are  avoided.  Moreover,  the  component  modes  will  experience  the  same 
tuning  due  to  bound  electrons,  so  that  perturbations  due  to  atoms  and  ions 
cancel  and  the  method  should  work  well  in  weakly  ionized  plasmas.  Another 
feature  is  that,  unlike  other  methods,  there  is  in  principle  no  need  to  switch 
off  the  plasma  to  establish  a  reference  phase  or  frequency,  so  that  the  method 
is  appropriate  for  steady-state  plasmas,  for  example  MHD  flows,  which  cannot 
easily  be  switched  on  and  off.  On  the  other  hand  the  low  beat  frequency  would 
preclude  fast  time  resolution  with  pulsed  plasmas.  However,  this  may  be  over¬ 
come  by  using  a  suitable  frequency  off-set  or  bias  as  described  below. 

It  is  clear  that  to  allow  the  linearly  polarized  mode  to  rotate,  the 
laser  must  have  rotational  symmetry.  Thus  Brewster  windows  are  precluded;  all 
windows  must  be  normal  to  the  axis  and  should  be  anti-reflection  coated.  It 
also  appears  necessary  to  ensure  operation  in  a  single  axial  mode.  While  dif¬ 
ferent  axial  modes  will  give  the  same  splitting  6f,  so  that  the  composite  modes 
rotate  at  the  same  frequency,  their  planes  of  polarizations  may  be  different 
and  give  rise  to  detected  signals  at  multiples  of  6f. 

It  is  anticipated  that  the  minimum  measurable  frequency  splitting  6fm^n  will 
be  set  by  the  tendency  for  closely  spaced  modes  to  lock  together  as  a  result  of 
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stray  coupling.  The  situation  is  similar  to  that  in  ring  lasers  [10],  where 
scattering  at  the  mirrors  causes  modes  propagating  in  opposite  directions  around 
the  ring  to  lock  for  roations  giving  frequency  shifts  less  than  100-1000  Hz.  In 
the  present  case  scattering  at  imperfections  of  the  mirrors  or  windows  may  be 
expected  to  couple  the  component  modes  and  lock  the  composite  mode's  plane  of 
polarization  onto  the  imperfections.  Judging  by  ring  laser  experience  this  may 
be  estimated  to  occur  for  splittings  less  than  -1  kHz. 

Mode  locking  may  be  prevented  by  producing  a  frequency  bias,  either  by  intro¬ 
ducing  an  additional  splitting  with  a  Faraday  rotator,  or  by  Zeeman  splitting  of 
the  laser  modes  [11]  on  application  of  a  magnetic  field.  By  biasing  to  a  suitably 
high  frequency  and  using  a  fast  FM  discriminator  it  appears  that  temporal  reso¬ 
lution  of  density  changes  may  be  obtained  for  pulsed  plasmas. 

The  feasibility  of  the  proposed  method  depends,  on  the  one  hand,  on  the 

configuration  and  parameters  f^,  f£  of  the  plasma,  and,  on  the  other,  on  the 

availability  of  suitable  lasers.  Since  the  optical  axis  must  lie  along  the 

magnetic  field  the  method  is  suited  for  systems  with  axial  symmetry  such  as 

mirror  machines  and  for  MHD  flows.  For  example,  in  the  combustion  MHD  facility 
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at  Stanford  the  duct  width  £  *  10  cm,  B  s  25  kG  and  n^  s  10  cm  .  Considering 
a  He-Ne  laser  of  mirror  spacing  L  =  100  cm,  the  frequency  shifts  6f  calculated 
from  Eq.  6  for  the  three  commonly  used  Ne  transitions  are  tabulated  below. 


x  M 

ftf  (kHz) 

0.633 

0.25 

1.152 

0.82 

3.39 

T.i 

For  these  conditions  the  0.633  ym  line  gives  a  rather  small  shift  and  in  any 
case  might  give  difficulties  because  of  its  low  gain  in  view  of  irreducible 
reflection  losses  at  the  windows.  However,  either  the  1.152  ym  line  or  the  high 
gain  3.39  ym  line  appear  to  be  convenient. 

Finally  it  may  be  noted  that  for  a  plasma  magnetized  perpendicular  to  the 
laser  axis  one  should  get  splitting  into  orthogonal  linearly  polarized  modes 
corresponding  to  the  ordinary  and  extraordinary  plasma  waves.  However,  in  that 
case,  the  splitting  is  smaller  again  by  a  factor  (f c/ f ) . 
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4.2  Velocity  measurement  by  a  laser  doppler  velocimeter 


4.2.1  Introduction 

Measurement  of  the  doppler  shift  of  focused  laser  light  offers  a 
powerful  non- perturbing  technique  for  probing  various  types  of  fluid  flows. 
Recently  it  has  been  used  in  a  variety  of  flows  including  low-velocity 
liquid  flows,  hypersonic  wind  tunnels,  rocket  exhausts,  and  also  in  high 
temperature  highly- ionized  plasmas  for  thermonuclear  fusion  research.  In 
principle  it  can  yield  information  with  high  spatial  resolution,  on  the 
average  velocity  and  the  turbulent  velocity  spectrum  o:  the  scattering 
centers.  In  the  case  of  neutral  or  weakly- ionized  gas  flows,  the  effec¬ 
tive  scattering  centers  are  macroscopic  particles  in  the  size  range  1-10 
micron  which  are  either  naturally  present  in  the  flow  or  deliberately  added 
to  facilitate  the  measurement. 

For  MHD  flows,  great  interest  attaches  to  the  structure  of  the 
gas-dynamic  boundary  layer  at  the  electrodes,  and  its  modification  by  MHD 
action  due  to  collisional  coupling  between  the  neutral  gas  and  the  plasma. 
The  study  of  this  boundary  layer  presents  a  considerably  more  challenging 
problem  than  appears  to  have  been  met  hitrrto  using  laser  doppler  tech¬ 
niques  in  other  types  of  flows.  The  requirement  on  spatial  resolution  is 
particularly  stringent.  For  typical  conditions  in  the  Stanford  combustion 
M.H.D.  channel,  the  boundary  layer  thickness  is  yQ  ~  1  cm  ,  while,  on  the 
assumption  of  a  seventh  power  velocity  profile  (u/uq)  =  (y/yQ)^^  .  the 
position  where  the  velocity  u  is  half  the  free-stream  value  uq  ,  is 
only  y  —  0.1  mm  (E  100  y)  from  the  surface.  Thus,  to  obtain  useful  data 
on  the  boundary  layer,  the  spatial  resolution  should  be  comparable  to  or 


22 


better  than  this,  and  the  instrument  must  be  capable  of  measuring  the  flow 
velocity  to  within  this  distance  of  an  incandescent  electrode.  Other  prob¬ 
lems  involve  obtaining  optical  access  to  the  measurement  point  without 
disturbing  the  flow  or  exposing  the  window  to  damage  from  the  hot  gas 
stream,  ensuring  an  adequate  concentration  of  scattering  centers  of  suitable 
size  and  elimination  of  the  effects  of  vibration  and  strong  magnetic  fields. 

In  the  following  sections  we  report  the  design  and  preliminary  testing 
of  a  laser  doppler  velocimeter  to  satisfy  these  requirements. 

4.2.  2.  Basic  Design  Considerations 

In  this  section  we  outline  the  basic  considerations  leading  to  the 
choice  of  configuration  discussed  in  detail  later. 

(a)  Scattering  Geometry 

The  first  basic  choice  concerns  the  scattering  geometry. 


Particle  Scattering  Geometry 
Figure  8 

For  the  general  geometry  of  Fig.  8,  the  frequency  shift  5  v  is  given 
by 

^  =  2(f)  cos  e  sin  |  .  (1) 
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For  Instance,  for  visible  light  v  ~  6  x  lO*-^  Hz  ,  and  taking  cos  9  sin(£/2) 
~  l  ,  then  for  a  velocity  v  ~  500  m/sec  ,  typical  of  the  Stanford  combus¬ 
tion  MHD  rig,  we  have  5v  ~  2  GHz,  a  microwave  frequency. 

In  practice,  the  backscatter  configuration,  ^  =  rr  ,  is  very  convenient, 
indeed  virtually  mandatory,  since  the  same  optics  are  used  for  the  incident 
and  scattered  beams.  In  this  way,  alignment  of  these  beams  at  a  common 
point  is  automatically  ensured,  thus  avoiding  the  difficult  tracking  prob¬ 
lem  involved  in  using  separate  optics  for  the  two  beams.  The  backscatter 
configuration  is  also  economical  on  optical  components  and  requires  only  a 
single  viewing  port  into  the  MHD  channel.  Furthermore,  the  distribution 
of  scattered  power,  according  to  the  theory  of  Hie  scattering  from  spheres, 
shows  a  relative  maximum  in  the  backscatter  direction.  On  the  debit  side 
is  the  fact  that  the  incident  and  scattered  beams  must  be  separated  by  a 
beam  splitter,  which  wastes  half  the  power  in  each  beam.  However ,  this 
power  loss  is  acceptable.  Another  problem  encountered  with  the  backscat¬ 
ter  configuration  is  that  the  power  incidentally  scattered  and  reflected 
back  from  the  incident  beam  by  the  optical  components,  may  give  rise  to 
a  reference  signal  at  the  unshifted  frequency  which  is  too  large  compared 
with  the  signal  scattered  from  moving  particles.  This  difficulty  may  be 
overcome  by  polarization  decoupling,  as  described  later. 

(b)  Frequency  Shift  Measurement 

There  are  two  basic  methods  for  measuring  the  frequency  shift:  optical 
heterodyning  and  direct  spectral  analysis  using  a  high  resolution  Fabry- 
Perot  interferometer.  The  former  is  relatively  insensitive  tc  frequency 
jitter  and  drift  due  to  vibration  and  temperature  changes,  since  the  scat¬ 
tered  signal  is  instantaneously  mixed  with  a  fraction  of  the  incident  signal. 
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In  contrast  the  second  nethod  requires  a  high  degree  of  frequency  stability 
of  both  the  laser  and  the  inter fer one ter.  Unfortunately,  optical  hetero¬ 
dyning  is  limited  by  the  frequency  response  of  photomultipliers  to  shifts 
less  than  500  MHz,  thus  restricting  its  use  to  low  velocities.  One  pos¬ 
sible  way  of  circumventing  this,  in  the  backscatter  mode,  is  to  look  almost 
normal  to  the  flow,  i.e.,  8  ~  90®.  However  this  is  undesirable  because  the 
system  is  then  relatively  mote  sensitive  to  turbulent  velocity  components 
along  the  line  of  sight.  Furthermore,  with  finite  apertures,  the  range  of 
angles  about  9  90°  causes  serious  degradation  of  velocity  resolution. 

Possibly  optical  heterodyning  to  frequencies  as  high  as  10  GHz  could  be 
accomplished  using  the  latest,  state-of-the-art  fast  photo- detectors.  How¬ 
ever  these  have  low  gain  and  would  require  wide-band  amplifiers.  For  these 
reasons,  it  was  decided  to  employ  the  second,  direct  method,  and  assume 
that  the  necessary  frequency  stability  can  be  achieved. 

The  conventional  Fabry-Perot  interferometer  consists  of  a  pair  of 
air-spaced,  highly- re flee ting,  optical  flats,  aligned  parallel  to  each 
other  and  normal  to  the  light  beam.  It  acts  as  a  narrow  band  frequency 
filter,  transmitting  collimated  light  only  when  the  plate  separation  t 
and  wavelength  \  satisfy  the  condition  for  constructive  interference  of 
the  multiply  reflected  beams: 

2t  =  n\  .  (2) 

The  frequency  interval  between  successive  modes,  n  ,  n+1  ,  is  called  the 
free  spectral  range  Av  given  by: 

Av  =  c/2t  .  (3) 
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This  is  the  maximum  frequency  difference  that  can  be  measured  without  am¬ 
biguity  due  to  the  multimode  character  of  the  device.  The  width  d-„  to 
half- power  of  each  mode  is  given  by 


dv  =  £v/:t  ,  (A) 

where  Y  is  the  finesse  and  is  determined  by  the  flatness  of  the  plates 
and  the  reflectivity  of  the  coatings.  With  high  quality  optical  flats  and 
coatings, values  of  •*  -  100  are  attainable. 

For  the  present  purpose,  we  would  choose  the  spacing  t  to  give  a 
free  spectral  range  somewhat  greater  than  the  maximum  frequency  shift  cor¬ 
responding  to  the  maximum  velocity  to  be  measured.  For  instance  t  7.5  cm 
gives  Av  2  GHz  .  By  arranging  for  a  reference  beam,  as  well  as  the 
scattered  beam,  to  enter  the  interferometer,  and  by  scanning  the  plate 
separation  over  a  distance  somewhat  more  than  half  a  wavelength,  transmitted 
signals  may  be  displayed  on  a  synchronously  driven  oscilloscope  as  shown  in 
Fig.  9.; 
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FREQUENCY  OR 
VELOCITY  SCALE 


Figure  9 

Oscilloscope  Velocity  Display 
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The  horizontal  deflection  corresponds  to  a  linear  frequency  scale  for  the 
frequency  shift  *v  of  the  scattered  light,  with  reference  Barkers  at 
zero  frequency  and  at  \v  ,  the  latter  being  calculable  from  the  known 
plate  separation.  It  should  be  noted  that  if  -v  exceeds  £y  ,  the 
signal  will  appear  at  a  position  corresponding  to  the  frequency  (Viv)  , 
sc  that  r'.iere  is  ambiguity  in  the  freqeucny  shift  to  within  an  integral 
number  ol  £v  -  This  is  the  reason  why  it  is  desirable  to  choose  Av 
somewhat  larger  than  the  maximum  expected  frequency  shift.  The  velocity 
resolution,  expressed  as  a  percentage  of  the  maximum  velocity  correspond¬ 
ing  to  »  is  then  (100/.t)7.  ,  so  that  a  velocity  resolution  of  the  order 
of  n.  is  possible. 

Recently,  piezo-electr'.cally  scanned  Fabry-Perot  interferometers 
using  spherical  mirrors  spaced  confocally,  have  become  conmercially  avail¬ 
able.  These  have  a  number  of  advantages  over  the  planar  type,  including 
lower  cost,  higher  finesse  and,  most  importantly,  a  greater  light  gathering 
power  (dtendue)  when  used  wich  well- collimated  narrow  beams.  This  latter 
feature  is  particularly  appropriate  for  the  optical  configuration  considered 
below,  in  which  the  scattered  beam  is  compressed  by  reverse  jassage  through 
the  same  telescope  that  is  used  to  expand  the  incident  beam.  The  only  dis¬ 
advantage  of  thf  confocal  Fabry-Perot  is  that  the  free- spectral  range  is 
not  adjustable,  being  fixed  for  a  given  pair  of  mirrors  by  the  confocal 
spacing  condition.  However,  various  free  spectral  ranges  in  the:  range 
1-10  GHz  can  be  obtained  by  inserting  different  mirror  sets  into  the  same 
housing.  Typically  a  finesse  of  300  is  achieved  in  the  visible,  corres¬ 
ponding  to  a  velocity  resolution  of  0.37«  of  the  maximum  velocity.  This  is 
more  than  adequate  for  our  purpose. 
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In  the  case  of  the  confocal  Fabry-Perot,  formulae  (2)  and  (3)  are 


modified  to 


4  r  =  n). 

(2a) 

and 

$ 

ii 

n 

<r 

(2b) 

where  r 

is  the  radius  of  curvature  of  the  mirrors  and  equals  their 

spacing. 

(c) 

Beam  Focusing  and  Spatial  Resolution 

Lasers  operating  single-mode  in  the  visible  spectrum  typically  produce 

2 

an  output  beam  with  a  gaussian  intensity  profile  of  width  to  1/e  points 

—  1  ran  and  with  spherical  wavefronts  corresponding  to  a  divergence  half- 
_  3 

angle  of  —10  radian.  To  obtain  the  required  spatial  resolution  this 
beam  must  be  focused. 

Now  if  a  gaussian  beam  is  focused  by  an  aberration- free  lens  (or 
mirror),  then  the  intensity  profile  remains  gaussian  and  is  given  by: 


where 

by: 


P 

o 


2  P 

I(r,z)  = - s -  exp  -  2 

tt  aZ(z) 

is  the  total  power  and  the  radius 


2 

a  (to  l/e“  point)  is  given 


,  N2_l/2 

»<*>  -  »,[i  +  r*-?)  ]  •  <« 

L  \  a  /  -i 

s 

Here  z  is  measured  from  the  focal  plane  where  the  beam  has  its  minimum 
radius  a 

s 

If  a  collimated  beam  of  dimeter  2  a^  is  focused  by  a  lens  of  focal 
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length  i  ,  as  shown  In  Fig .  10 


we  have  from  Eq.  (6): 


Figure  10 

Beam  Focusing  Diagram 


and 
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(8) 


where  £  is  the  axial  distance  to  the  half- intensity  points  and  F  = 
s 

f/2ab  is  the  effective  numerical  aperture  of  the  lens.  It  should  be  noted 
that  away  from  the  focal  spot,  the  intensity  falls  off  much  less  rapidly 
in  the  axial  direction  than  in  the  radial  direction,  and,  except  in  the 
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limit  of  F/l  optics,  the  length  of  the  focal  region  is  greater  than  its 
diameter  (by  a  factor  — 2F  ). 

Physical  restrictions  set  by  the  geometry  of  the  Stanford  M.H.D.  rig 
prevent  our  putting  a  focusing  lens  closer  than  about  50  cm  from  the  obser¬ 
vation  point,  i.e.  f  >  50  cm  .  If  the  laser  beam  were  directly  focused 
by  a  lens  of  this  focal  length,  we  should  have  F  ~  500  ,  and  hence  a  focal 

region  of  diameter  2  a  —  0.3  mm  and  length  2ft  ~  15  cm  .  This  is  much 

s  s 

too  large  for  our  purpose  and  we  are  led  to  consider  expanding  the  beam  to 
a  larger  diameter  before  focusing  it.  Diffraction- limited  (i.e.  aberration 
free)  beam-expanding  telescopes  for  use  with  visible  lasers  are  commer¬ 
cially  available  with  expansion  factors,  depending  on  choice  of  lenses,  up 
to  50  times.  They  consist  of  an  input  beam- expanding  lens  (effectively  a 
microscope  objective)  and  a  large  diameter  collimating  lens  combination 
or  telescope  objective  designed  for  minimum  spherical  observation. 

Using  such  a  beam- expanding  telescope,  our  basic  optical  configuration 
then  appears  as  shown  in  Fig.  11. 


Figure  H 

Optical  Configuration 
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It  should  be  noted  that  the  focusing  lens  also  acts  as  the  collecting  lens 
for  scattered  light,  which  is  compressed  into  a  narrow  beam,  —  1  mm  diame¬ 
ter,  for  acceptance  by  the  Fabry- Perot.  The  choice  of  beam  expansion  fac¬ 
tor,  and  hence  the  numerical  aperture  of  the  focusing  lens,  determines 
not  only  the  dimensions  of  the  focal  volume,  and  hence  the  spatial  resolu¬ 
tion,  but  also  the  collection  efficiency.  The  latter  improves  as  the  aper¬ 
ture  is  increased  (  F  decreased),  but  a  limit  is  set  by  the  degradation  in 
velocity  resolution  produced  by  the  range  of  angles  0  involved  in  the 
finite  aperture. 

Differentiating  Eq.  (1)  gives  the  relative  spread  in  apparent 
velocity,  due  to  the  finite  aperture  as: 

=  tan  g  d0  =  lanj  <  (9) 

The  angle  0  between  the  beam  axis  and  the  channel  axis  is  determined  by 
geometrical  considerations  as  0  =  26°  ,  as  discussed  later. 

In  practice,  the  beam  expansion  factor  was  chosen  to  be  a  factor  -50  , 

the  maximum  readily  available,  corresponding  to  an  expanded  beam  diameter 
of  50  mm  and  a  numerical  aperture  of  F  10  for  the  focusing  lens.  This 
is  consonant  with  considerations  of  the  geometry  of  the  window  and  port 
into  the  M.H.D.  channel.  Eq.  (6)  then  leads  to  an  apparent  frequency 
spread,  due  to  finite  aperture,  of  —  2.5%  .  In  practice,  taking  account 
of  the  energy  integrated  over  the  focusing  lens  area,  we  should  find  a 
relative  velocity  resolution  somewhat  less  than  this,  say  2 7»  at  most.; 

The  dimensions  of  the  focal  volume  calculated  from  Eqs.  (7)  and  (8) 
are  then  2a  7[j  ,  -  140jj  ,  which  should  give  excellent  spatial  rcsolu- 

tion  for  our  purposes.  It  may  be  noted  that  the  velocity  resolution  could 


31 


be  improved,  at  the  expense  of  u  duced  spatial  resolution,  by  using  a 
smaller  beam  expansion  factor  (by  changing  the  expanding  lens  in  the  tele¬ 
scope)  ,  but  it  should  be  remembered  that  the  signal  pulse  amplitude,  its 
duration  and  the  rate  of  acquisition  of  scattered  signals  are  all  affected 
by  such  changes . 

The  final  matter  to  be  discussed  under  the  heading  of  spatial  resolu- 
tic  is  spatial  filtering.  By  placing  a  pinhole  at  the  focus  of  the  beam¬ 
expanding  lens,  i.e.  at  the  cross-over  in  the  telescope,  one  can  discrimi¬ 
nate  against  the  detection  of  light  originating  outside  the  focal  volume  of 
the  focusing  lens.  This  not  only  sharpens  the  ;patial  resolution  for  scat¬ 
tered  laser  light  (especially  in  the  axial  direction),  but  discriminates 
against  the  detection  of  radiation  from  the  hot  gas  and,  more  especially, 
from  the  hot  electrode,  particularly  when  probing  very  close  to  the  latter. 
Rejection  of  the  radiation  background  is  also  achieved  by  using  a  narrow 
band  filter  (~  10 A. U. centered  on  the  laser  line)  in  front  of  the  detector, 
but  spatial  filtering  is  also  effective  and  desirable.  Spatial  filters, 
in  the  form  of  a  pinhole  with  accurate  x  ,  y  ,  z  position  controls,  are 
supplied  as  accessories  to  beam- expanding  telescopes.  By  using  a  pinhole 
somewhat  larger  than  the  theoretical  cross-over  diameter,  one  can  clean  up 
any  wavefront  distortion,  either  inherent  in  the  laser  output,  or  intro¬ 
duced  by  imperfections  in  the  beam- expanding  lens.  Thus  the  spatial  filter 
serves  the  double  purpose  of  maintaining  wave- front  purity  in  the  incident 
beam,  and  discrimination  against  signals  from  outside  the  focal  volume  in 
the  collected  beam. 

(d)  Laser  Source  Requirements 

For  proper  functioning  of  the  velocity  analyzer,  it  is  clearly  neces¬ 
sary  that  the  laser  source  should  operate  single-mode.  Furthermore  it  is 
clear  that  because  of  the  very  small  focal  volume,  necessitated  by  the  high 
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spatial  resolution  required,  then,  for  any  realizable  particle  concentration 
in  the  flow,  the  probability  qf  finding  a  particle  in  the  focal  volume  at 
any  instant  is  much  less  than  unity.  Consequently  we  have  to  detect  indi¬ 
vidual  pulses  scattered  by  single  particles  as  they  pass  through  the  focal 

volume.  For  instance,  for  a  focal  volume  of  cross-sectional  area  projected 

-  5  2 

normal  to  the  flow  of  10  cm  ,  a  flow  velocity  of  500  m/sec,  and  a  particle 

3-3 

concentration  of  10  cm  (1  mm  average  spacing  between  particles) ,  we 
shall  have  *“  500  particles  passing  through  the  focal  volume  per  second. 
However,  because  the  Fabry-Perot  scans  with  a  frequency  (velocity)  band¬ 
width  only  ~  17»  of  the  frequency  (velocity)  excursion,  we  shall  detect  only 
~  17«  of  these,  i.e.  ,  *“  5  particles/sec  for  these  conditions.  , 

•  i 

The  problem  thus  comes  down  to  the  detectability  of  such  pulses  against 
the  noise  background  which  in  this  case,  assuming  negligible  background 
radiation  is  admitted  to  the  detector,  consists  of  the  dark  current  pulses 
of  the  photomultiplier.  The  amplitude  of  the  signal  pulses  is  proportional 
to  the  laser  power  level  and  the  scatteringcross-section  of  a  particle,  as 
well  as  the  collection  and  transmission  efficiencies  of  the  optics.  Clearly, 
in  view  of  the  low  rate  of  acquisition  of  signals,  it  is  desirable  to  use  a 
continuously  operating  laser,  rather  than  a  pulsed  one,  provided  sufficient 
power  can  be  obtained  in  the  continuous  mode  to  give  detectable  signals 
from  the  size  of  particles  considered. 

Since  tne  quantitative  evaluation  of  these  factors,  involved  in  signal 
detectability,  is  subject  to  considerable  uncertainty,  it  is  best  to  rely  on 

I 

performance  data  from  previous  work.  With  a  laser  doppler  velocimeter 
rather  similar  to  that  considered  here,  Seifert  was  able  to  detect  signals 
from  particles  in  the  size  range  1  -  10  (a  in  solid  fuel  rocket  exhausts, 
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using  a  continuous,  single-mode  He-Ne  laser  with  an  output  of  ~100  |i  watts. 
However,  at  the  lower  size  limit  (~1  p)  the  signal/noise  ratio  was  marginal 
and  he  had  to  resort  to  dividing  the  scattered  light  between  two  photomul¬ 
tipliers  and  using  a  coincidence  circuit  to  improve  the  detection  of  sig¬ 
nals  against  the  random  dark  current  pulse  background. 

For  these  reasons  it  seemed  desirable  to  use  an  Argon  ion  laser,  which 
is  the  only  convenient  alternative,  and  which  delivers  some  100  m  watts 

single  mode  on  either  a  green  line  (5145  Au)  or  a  blue  line  (4880  Au) . 

3 

This  gives  a  factor  10  in  power  compared  with  the  He-Ne  laser  and  further¬ 
more,  one  also  gains  from  the  higher  quantum  efficiency  of  photocathodes  in 
the  blue-green  compared  with  the  red. 

The  only  disadvantage  of  the  A  laser  compared  with  the  He-Ne  one, 
is  that,  because  of  the  much  higher  discharge  power,  the  amplitude  and 
frequency  stability  is  somewhat  poorer.  The  principal  effect  is  low  fre¬ 
quency  (a  few  Hz)  modulation  of  the  laser  frequency  due  to  mirror  vibra¬ 
tion  induced  by  turbulence  or  local  boiling  in  the  water  cooling  jacket. 

The  frequency  jitter  is  of  the  order  +  20  MHz,  which  causes  some  degrada¬ 
tion  of  the  velocity  resolution  due  to  jitter  (along  the  horizontal  axis) 
in  the  display  (Fig.  2)  over  times  of  the  order  of  0.1  secs.  Whether  this 
is  important  depends  on  how  long  the  display  must  be  stored  before  record¬ 
ing,  which  in  turn  depends  on  the  signal  acquisition  rate.  If  the  latter 
is  high  enough,  the  display  can  be  stored  for  a  time  short  compared  with 
the  characteristic  jitter  time,  in  which  case  the  jitter  is  unimportant. 

(e)  Detection,  Display  and  Recording 

The  signals  to  be  detected  consist  of  reference  signals  and  scattered 
signals  whose  characteristics  we  consider  in  turn.  Ideally,  the  amplitude 
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of  the  reference  signals  should  be  of  the  same  order  or  somewhat  larger 
than  the  scattered  signals,  which  themselves  should  be  at  least  an  order  of 
magnitude  than  the  noise  level.  On  first  testing  the  system  in  the  back- 
scatter  configuration  shown  in  Fig.  H  it  was  found  that  the  reference 
signal  due  to  reflection  at  the  (anti-reflection  coated)  surfaces  of  tie 
beam- expanding  lens,  together  with  light  scattered  from  the  edges  of  the 
pinhole,  constituted  a  reference  signal  which  was  much  too  large — some 
three  orders  larger  than  noise.  To  reduce  this  level,  without  affecting 
the  scattered  signal  level,  use  was  made  of  a  polarization  decoupling  scheme. 
A  horizontally-oriented  polarizer  was  introduced  in  front  of  the  Fabry- 
Perot.  Since  the  laser  output  is  vertically  polarized  this  rejects  the 
reference  signal.  A  quarter  wave  plate,  designed  for  \  =  5145  Au  ,  was 
placed  between  the  telescope  objective  and  the  focusing  lens,  so  oriented 
to  produce  circularly  polarized  light  in  the  probing  beam.  Scattered 
light  which  is  likewise  circularly  polarized,  is  converted  to  a  horizon¬ 
tally  polarized  signal  on  reverse  passage  through  the  quarter  wave  plate, 
and  passes  the  horizontal  polarizer  to  be  detected.  In  this  way  the  ampli¬ 
tude  of  the  reference  beam  was  reduced  by  some  20  db  (i.e.  two  orders) 
relative  to  the  signal. 

The  pulse  length  of  the  reference  signal  is  determined  by  the  band¬ 
width  of  the  Fabry-Perot  together  with  its  scanning  rate.  The  scan-rate 
cannot  exceed  about  1  kHz  if  the  mirrors  are  to  faithfully  follow  the 
sawtooth  waveform,  and  in  practice  a  time  of  ^-10  m  sec  to  scan  over  one 
free-spectral  range  was  found  to  be  convenient..  For  a  finesse  of  100, 
this  gives  a  reference  signal  pulse  length  of  -100  ^  sec.  Such  pulses  will 
write  quite  satisfactorily  on  a  synchronously  driven  scope  with  a  time  base 
of  10  m  sec 
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By  comparison,  the  scattered  light  pulse  due  to  passage  of  a  single 
particle  through  the  focal  volume  is  of  very  short  duration,  —  50  n  sec 
for  the  conditions  envisaged,  and  at  the  lower  limit  of  particle  size 
(<  1  u)  may  contain  only  of  the  order  of  10  photons.  The  most  convenient 
and  versatile  detector  is  a  photomultiplier,  since  it  is  fast  and  has 
high  gain,  capable  of  detecting  such  pulses  and  amplifying  them  to  a  level 
at  which  they  can  be  recorded. 

For  instance,  the  14- stage  RCA  7205  with  S.20  cathode  has  a  quantum 
efficiency  of  12%  at  5145  Au,  a  rise  time  of  3  n  sec,  and  a  current  gain  of 
5  x  10^  when  operated  at  2400  V.  Thus,  using  a  50  £2  output  cable  terminated 
in  50  ,  one  can  expect  to  reproduce  pulses  as  short  as  10  n  sec,  con¬ 

sisting  of  only  a  few  photoelectrons,  and  amplify  them  to  a  level  —100  m  V, 
which  is  adequate  as  input  to  a  broad-band  (30  MHz)  oscilloscope  amplifier. 
It  is  desirable  to  be  able  to  record  individual  scattered  light  pulses 
without  loss  of  time  resolution,  since  it  gives  an  independent  check  on 
the  effective  size  of  the  focal  volume  actually  achieved.  However  it  is 
generally  difficult  or  inconvenient  to  record  such  fast  pulses  photo¬ 
graphically  from  an  oscilloscope  and,  except  when  investigating  the  pulse 
length,  it  is  desirable  to  lengthen  the  pulses  so  that  they  write  more  ef¬ 
fectively  on  the  oscilloscope  screen.  This  can  be  done  mort  easily  by 
using  a  much  larger  load  resistor,  so  that  the  time  constant  formed  with 
the  output  capacity  of  the  photomultiplier,  together  with  that  of  the  cable, 
is  long  enough  for  effective  writing.  At  the  s  ime  time  the  photomultiplier 
gain  is  reduced  by  lowering  the  operating  voltage,  to  avoid  saturation 
ef  fee  ts  .; 

For  a  velocity  display  as  shown  in  Fig.-  9,  it  is  desirable  to  increase 


the  detector  time  constant  until  it  is  comparable  with  the  pulse  length 
(—100  |j  sec)  cf  the  reference  signals  as  given  by  the  Fabry- Perot  band¬ 
width  and  scan  refe.  The  writing  efficiency  is  then  comparable  for  the 
scattered  and  reference  signals  when  these  are  of  comparable  amplitudes. 

The  noise  background,  against  which  the  s'gnal  pulses  must  be  detected, 
consists  of  the  so-called  random  dark  currant  pulses.  These  arise  from 
thermonic  emission  from  the  cathode,  from  secondary  emission  from  the  elec¬ 
trodes  due  to  ion  bombardment,  and  from  cold  emission  from  electrodes. 

While  the  pulses  originate  primarily  from  single  electron  emission  events, 
the  anode  pulses  vary  in  height  due  to  the  varying  statistics  of  the  multi¬ 
plication  process.  The  average  frequency  of  such  pulses  is  not  normally 
quoted,  but  may  be  estimated  from  the  manufacturer's  specification.  From 

o 

the  average  anode  dark  current  (at  22  C)  and  the  gain,  we  obtain  an  equiva- 

-14  4 

lent  cathode  dark  current  —10  A  ,  corresponding  to  some  5  x  10  elec¬ 
trons/sec  randomly  generated  at  the  cathode.  The  frequency  of  such  pulses 
may  be  reduced  by  cooling  the  photomultiplier  with  dry  ice,  which  elimi¬ 
nates  the  thermonically  emitted  component,  but  there  is  an  irreducible 

4 

minimum  pulse  rate  of  —10  /sec  . 

To  detect  signal  pulses  with  reasonable  certainty  against  this  back¬ 
ground,  either  they  must  be  significantly  larger  than  the  dark  current 
pulses  or  they  must  occur  significantly  more  frequently.  We  have  seen 
that,  in  view  of  the  small  focal  volume,  it  is  unlikely,  for  reasonable 
particle  concentrations,  that  we  shall  have  a  signal  rate  of  more  than 
—10  per  second.  Consequently  the  second  alternative  is  not  open  to  us, 
and  it  is  essential  that  the  signal  pulses  be  larger  than  the  dark  current- 
pulses,  containing,  say,  at  least  100  photons  (12  photo  electrons)  per 
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pulse.  This  sets  a  limit  on  how  small  a  particle  can  be  detected. 

In  section  4  we  make  some  rough  calculations  of  the  amplitude  of  the 

detected  signals  for  our  particular  instrument  design,  and  show  that  we 
3 

may  expect  -JLO  photons  per  pulse  for  1  p  particles.  It  should  be  empha¬ 
sized,  however,  that  this  is  a  very  rough  estimate  and  that  the  performance 
of  the  system  can  only  be  reliably  determined  from  measurements  using 
particles  of  known  size. 

If,  for  some  reason,  the  signals  are  not  so  much  larger  than  the  dark 
current  pulses  that  they  can  be  unambiguously  distinguished,  then  it  is 
still  possible  to  distinguish  them,  even  though  they  occur  much  less  fre¬ 
quently  than  the  dark  current  pulses,  by  dividing  the  signal  between  two 
photomultipliers  with  fast  response,  followed  by  a  coincidence  circuit. 
Hopefully  our  signals  will  be  large  enough  that  we  shall  not  need  to  re¬ 
sort  to  this  degree  of  sophistication.  The  primary  mode  of  operation  uses 
a  velocity  display  as  shown  in  Fig.  9.  Since  the  scan  period  is  -10  m  sec, 
while  the  signal  rate  is  anticipated  to  be  only  ~10  per  second,  it  is  clear 
that,  in  general,  we  shall  need  to  scan  continuously  and  store  signals  for 
the  order  of  one  or  more  seconds.  This  may  be  done  photographically, 
using  a  suitably  long  exposure,  but  a  more  convenient  method  is  to  use  a 
storage  oscilloscope.  In  this  way  a  suitable  result  is  obtained  before 
recording  it  permanently  on  film. 

A . 2 . 3 .  Detailed  Design  Considerations 

Having  determined  the  basic  design  parameters  of  the  velocimeter,  we 
now  discuss  some  detailed  matters  related  to  constructing  an  instrument 
capable  of  making  the  required  measurements  on  the  Stanford  combustion 
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MHD  rig.  The  latter  imposes  quite  severe  constraints  cnr.^eming  access 
to  the  measurement  region,  vibration  isolation  and  magnetic  shielding. 

First,  we  consider  the  question  of  access  to  the  measurement  region. 

(a)  Access  Port,  Scan  Mirror  and  Focusing  Lens 

The  main  MHD  channel  is  some  60  cm  long  and  situated  in  a  uniform 
magnetic  field  of  some  3  Tesla  provided  by  a  large  electromagnet.  The 
channel  is  constructed  of  double-walled,  stainless  steel  incorporating  a 
water  cooling  jacket,  lined  with  MgO  firebrick,  and  has  internal  dimensions 
10  cm  high  by  3  cm  wide.  It  is  required  to  make  measurements  adjacent  to 
an  electrode  (3  cm  wide  x  2  cm  long)  in  the  top  surface  of  the  duct  and 
positioned  near  the  downstream  end  of  the  main  channel,  where  the  boundary 
layer  is  fully  developed. 

The  access  port  must  transmit  the  probing  beam  which  consists  of  a 
converging  conical  pencil  of  half-angle  ~3° ,  corresponding  to  F/10  optics. 

It  must  also  accommodate  the  beam,  without  vignetting,  as  it  is  deflected 
vertically  to  move  the  focal  point  from  the  electrode  surface  to  some  2  cm 
below  it  into  the  free  stream. 

The  design  adopted,  shown  in  plan  in  Fig. 12  gives  access  through  a  side 
wall  downstream  of  the  measurement  point.  This  minimizes  perturbation  of 
the  boundary  layer  at  the  electrode  by  the  discontinuity  in  the  side  wall. 
The  angle  &  ,  in  the  horizontal  plane,  between  the  beam  and  channel  axes 
was  set  at  26°  by  geometrical  considerations.  This  gives  a  reduction 
factor  cos  fl^-0.90  in  Eq.  (1)  relating  the  frequency  shift  to  the 
velocity.  Considering  the  vertical  section  through  the  beam  and  port 
axis,  shown  in  Fig. 13,  it  is  clear  that  the  most  critical  factor  is  to 
ensure  that  the  convergent  pencil  can  be  focused  at  the  center  of  the 
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Figure  12  Plan  View  of  Port 


Fig. 13.  Vertical  Section  through  Port 


electrode  surface  without  vignetting  anywhere  on  the  port  _  hannel  walls. 
This  was  accomplished  by  making  the  angle  t  ,  between  the  horizontal  and 
the  axis  of  the  beam  (when  focused  on  the  electrode)  slightly  greater  than 
the  half-angle  of  the  convergent  cone.  The  upper  edge  of  the  cone  then  just 
clears  the  electrode  and  firebrick  surfaces. 

The  length  of  the  port  was  made  as  long  as  possible  (—12  in)  consistent 
with  not  obtruding  into  the  space  required  to  locate  the  scan  mirror,  so  as 
to  protect  the  window  as  far  as  possible  from  the  hot  gas  flow.  The  port 
walls  were  constructed  of  a  double  layer  of  stainless  steel  spaced  by  a 
water  cooling  channel,  and  the  whole  assembly  welded  to  the  side  cover  of  the 
main  channel.  The  window,  consisting  of  a  2"  dia  x  1/2"  thick  fused  silica 
optical  flat  (\/20) ,  was  sealed  to  the  port  with  an  O-ring,  and  provision 
was  made  to  direct  a  stream  of  nitrogen  onto  its  inner  surface.  The 
nitrogen  flow  was  designed  to  inhibit  the  entry  of  hot  gas  into  the  port. 
However  the  flow  must  be  maintained  at  a  moderate  level  lest  it  perturb 
the  boundary  layer  at  the  measurement  point. 

To  keep  the  diameters  of  the  scan  mirror  and  focusing  lens  of 
reasonable  size,  it  is  necessary  to  position  them  as  near  as  possible  to 
the  measurement  point.  Of  the  two  possible  configurations  of  scan  mirror 
and  focusing  lens,  corresponding  to  deflection  before  or  after  focusing, 
the  latter  arrangement  was  adopted,  since  the  former  would  involve  using  the 
lens  off-axis  and  thus  give  rise  to  coma  in  the  focal  "spot".  The  scan 
mirror,  consisting  of  a  3"  dia  x  3/4"  thick  glass  optical  flat  () / 20)  is 
positioned  and  oriented  so  as  to  deflect  the  beam,  incident  from  a  direc¬ 
tion  normal  to  the  channel  axis,  centrally  into  the  port.  Its  front  sur¬ 
face  is  dielectric-coated  for  maximum  reflectance  at  this  angle  at 
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>  =  5145  A.D.  The  poor  reflection  at  other  wavelengths  minimizes  the  col¬ 
lection  of  radiation  from  the  gas  or  walls.  Scanning  of  the  focal  spot  in 
the  vertical  direction  is  accomplished  by  tilting  the  mirror  mount  about  a 
horizontal  axis  in  the  plane  of  the  mirror  surface.  This  motion  is  con¬ 
trolled  by  a  micrometer  positioned  3  in  from  the  tilt  axis.  Since  the  dis¬ 
tance  from  the  mirror  to  the  focal  spot  is  —21  in,  there  is  a  lever  arm 
m-'  gnification  of  seven  times  between  the  movements  of  the  focal  spot  and  the 
-aicrometer. 

The  focal  length  of  the  focusing  lens  should  be  chosen  to  be  the  min¬ 
imum  consistent  with  locating  it  close  to  the  scan  mirror,  so  as  to  pre¬ 
serve  maximim  aperture.  This  resulted  in  the  value  f  =  56  cm,  which,  in 
view  of  the  maximum  diameter  (50  mm)  of  the  expanded  beam  set  by  the  tele¬ 
scope  aperture,  corresponds  to  F/ll  .  It  proved  impossible  to  find  a 
diffraction- limited  lens  (as  opposed  to  an  ordinary  achromat)  with  these 
approximate  characteristics,  except  as  a  special  order.  Eventually,  a 
60  mm  diameter  diffraction- limited  air-spaced  doublet  with  f  =  54  cm  (F/9) 
was  located.  This  focal  length  is  too  short,  and  with  a  collimated  beam 
from  the  telescope  would  put  the  focal  spot  some  2  cm  short  of  the  elec¬ 
trode  center.  However,  by  defocusing  the  telescope,  to  produce  a  slightly 
divergent  beam  (-0.1°  half -angle),  it  is  possible  to  bring  the  focal  spot 
to  the  required  positior  Whi1^  the  telescope  and  focusing  lens  are  de¬ 
signed  for  zero  spherical  abe  ration  for  collimated  beams,  the  degradation 
produced  by  this  slight  depa  ture  from  exact  collimation  is  quite  negligible. 

(b)  Layout  and  Mounts  g  of  Optical  System 

While  the  scan  mirror  and  focusing  lens  must  be  mounted  close  to  the 
port,  there  is  some  choice  possible  in  locating  the  rest  of  the  optical 
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system.  From  the  point  of  view  of  minimizing  the  magnetic  shielding  prob¬ 
lem,  and,  possibly,  also  the  vibration  isolation  problem,  it  would  be  de¬ 
sirable  to  mount  the  rest  of  the  system  on  a  stable  platform  well  removed 
(say  5  meters)  from  the  MHD  rig.  However,  the  operator  must  be  able  to 
set  and  read  the  scan  mirror  micrometer  as  well  as  make  adjustments  to  the 
laser  and  electronics  while  observing  the  oscilloscope  display.  Conse¬ 
quently,  a  remote  location  would  necessitate  an  accurate  remote  drive  and 
sensing  mechanism  for  the  scan  mirror.  For  these  reasons  a  compact  layout 
was  adopted,  as  shown  in  Fig.  14 ,  whereby  all  the  controls  are  within  reach 
of  a  single  operator. 

The  optical  system  is  mounted  on  a  rigid  framework  formed  of  two  4  in 
aluminum  H- beams,  bolted  at  right  angles  and  cantilevered  from  the  iron 
core  of  the  electromagnet.  Since  the  latter  weighs  some  10  tons  and  rests 
firmly  on  the  reinforced  concrete  floor,  it  was  judged  to  form  a  stable 

support.  The  free  ends  of  the  H- beams  are  supported  by  surveyor's  tripods. 
This  system  has  the  advantage  that  the  framework  can  be  readily  removed  and 
used  to  support  the  optical  system  in  an  alternative  location,  namely  a 
darkroom,  fo”  preliminary  testing.  With  the  exception  of  the  laser,  which 
rests  directly  on  the  H-beam,  all  the  optical  components,  including  the 
focusing  lens  and  scan  mirror,  are  mounted,  using  standard  optical  car¬ 
riages  and  adjustable  mounts,  on  an  optical  bench  which  is  rigidly  bolted 
to  the  H-beam.  This  greatly  facilitates  alignment  of  the  optics. 

It  was  necessary  to  screen  both  the  laser  and  the  photomultiplier 
from  the  stray  magnetic  field  from  the  electromagnet,  whose  magnitude  was 
-1000  oersted  at  each  location.  While  the  laser  operates  in  an  axial  field 
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Figure  14  Otical  Layout 


of  -800  oersted,  it  is  adversely  affected  by  transverse  fields  of  more  than 
a  few  oersted.  Adequate  shielding  was  provided  by  surroundii.,  the  laser 
with  an  open-ended  box  constructed  of  1/8"  thick  mild  steel.  The  function¬ 
ing  of  the  photomultiplier  is  likewise  adversely  affected  by  fields  of  a 
few  oersted.  Adequate  shielding  was  provided  by  surrounding  the  tube  by  a 
length  of  thick-walled  (1/4")  iron  pipe.  All  the  components  of  the  optical 
system  are  constructed  of  non-magnetic  material,  with  the  exception  of 
the  micrometers  controlling  the  tilt  motions  of  the  optical  bench  mounts. 

The  force  exerted  on  the  micrometers  could  possibly  cause  misalignment  of 
the  optics  during  operation  when  the  magnetic  field  is  changed.  In  practice, 
the  only  noticeable  effect  was  a  deflection  of  the  focal  spot  due  to  strain 
induced  in  the  mount  of  the  scan  mirror.  This  could  be  corrected  by  re¬ 
setting  the  miciometer  (by  0.004"),  but  a  better  solutic  would  be  to 
replace  the  micrometer  by  a  non-magnetic  one.. 
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4.2.4  Component  Specifications  and  System  Characteristics 


The  specifications  of  the  optical  components  selected  for  the  final 


design  are  listed  in  Table  I. 


TABLE  I 


Laser  Power,  Single  Mode  (5145  Au) 

P 

0 

* 

r*H 

I 

o 

H 

Output  beam  diameter 

2aL 

1.1  mm 

"  "  divergence 

■ 

9x10  f*  radian 

Transmission  of  Beam-Splitter 

tb 

507„ 

Telescope  expanding  lens ,  focal  length  , 

fE 

6.1  mm 

"  "  "  Aperture  ' 

2  ag 

1 . 7  mm 

Spatial  filter,  pinhole  diameter 

2aph 

15  U 

Telescope  collimating  lens,  focal  length 

fx 

200  mm 

"  "  "  aperturd 

2  aT 

50  mm 

Transmission  of  Telescope  and  Spatial  Filter 

T 

T 

70% 

Transmission  of  Quarter-wave  plate  1 

T 

Q 

85% 

Focusing  Lens,  focal  length 

fF 

54  cm 

"  "  aperture 

60  mm 

"  "  transmission  ' 

tf 

95% 

Scan  Mirror,  Reflectivity 

*M 

99% 

1 

Window  transmission 

Tw  1 

90% 

Numerical  Aperture  of  Collection  optics 

F 

11 

Fabry- Perot- -Free  Spectral  Range 

Av 

2  or  8  GHz 

"  "  Finesse 

,  .:f 

—300 

"  "  Resolution 

dv 

7  or  27  MHz 

"  "  Entrance  Pupil 

2  mm 

"  "  Scan  Period 

'T 

-2 

10  sec 

"  "  Transmission 

T 

FP 

50% 

Polarizer,  Max.  Transmission 

T 

P 

90% 

Filter  Transmission 

*, 

50% 

Photomultiplier,  Quantum  Efficiency 

'  T1 

12% 

"  Current  Gain  at  2400  V 

5  x  107 
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From  these  specifications,  together  with  some  nominal  values  for  the 
flow  conditions,  we  can  now  detail  some  calculations  of  the  anticipated 
system  characteristics.  The  figures  are  rough  estimates  and  probably 
optimistic,  particularly  with  regard  to  the  magnitude  of  the  detected 
scattered  signal.  However,  a  more  accurate  evaluation  is  hardly  feasible 
and  the  true  performance  must  be  determined  experimentally.  Nevertheless 
the  calculations  are  useful  to  illustrate  the  functional  dependence  of  the 
performance  on  the  external  parameters  of  the  design.  The  principal  re¬ 
sults  are  listed  in  Table  II.  Some  comments  on  the  various  entries  follow: 
(3,4,5)  The  particle  diameter  and  concentration  assumed  have  been  taken 

as  about  the  smallest  with  which  we  can  reasonably  work.  Under 

these  conditions,  the  volume  flow  of  particulate  matter  is  quite 

-  3  3  3 

low  (10  cm  /sec)  but  would  become  appreciable  (1  cm  /sec)  if 

6  -  3 

d  were  increased  to  10|i  cr  n  increased  to  10  cm  . 

(6)  The  input  beam  almost  exactly  fills  the  telescope  entrance  pupil. 

(7)  The  pinhole  was  chosen  to  be  five  times  the  theoretical  cross¬ 

over  diameter. 

(8,9)  The  telescope  objective  is  slightly  overfilled  and  the  expanded 
beam  diameter  is  determined  by  the  aperture  of  the  objective. 

(10)  By  ray  optics,  the  path  of  any  light  capable  of  being  detected 
(  must  pass  through  this  image  circle. 

(11,12)  The  focal  volume  is  represented  as  a  cylinder  with  uniform  il¬ 
lumination  within  it  and  zero  outside. 

(14)  The  angle  0  is  taken  when  the  focal  spot  is  at  the  electrode 
surface.; 

(17)  This  velocity  resolution  is  determined  by  the  finite  aperture  and 

not  by  the  bandwidth  of  the  Fabry-Perot,  which  is  relatively  smaller.; 


TABLE  II 


1. 

Channel  Cross-section 

A 

30  cm2 

2. 

Flow  Velocity 

V 

500  m/sec 

3. 

Particle  Diameter 

d 

lu  (10  ^*cm) 

4. 

Particle  Concentration 

n 

103cm'3 

5. 

Volume  Flow  of  Particulate 

Avnd3 

-3  3 

10  cm  /sec 

6. 

Beam  Diameter  at  Telescope 

/ 

Entrance 

2aL 

1.65  mm 

7. 

Diameter  of  Cross-over  in 
Telescope 

2  a 

X 

3^ 

8. 

Beam  Expansion  Factor 

M  =  fT^fE 

33 

9, 

Expanded  Beam  Diameter 

2%  =M.2aL 

50  mm 

10. 

Diameter  of  Pinhole  Image 

at  Focus 

2aph  “  2aph  fF/fT 

40m 

11. 

Diameter  of  Focal  Volume 

2as  =  (4/TT)fFx/2ab 

12. 

Length  of  Focal  Volume 

=  2tt  as2/X 

0. 15  mm 

13. 

Angle  between  beam  and  flow 
axes ,  horizontal 

e 

26’ 

14. 

Angle  between  beam  and  flow 
axes ,  vertical 

0 

3° 

15. 

Intercept  of  Focal  Volume 
with  flow  axis 

=  2ag/sin  a 

16u 

16. 

Spatial  Resolution  normal 
to  electrode 

fy  =  (2as+24  sin  0) 

15u 

17. 

Relative  Velocity  Resolution 

(*v/v)  =  tan  A/2F 

21 

18. 

Reference  Pulse  Length 

Tr  -  I  h 

30u  sec 

19. 

Scattered  Pulse  Length 

TS  =  6£/v 

30n  sec 

20. 

Area  of  Focal  «lume  normal 
to  flow 

A  =  2a  *2jg  sin  p 
s  s  s 

5  x  10  ^cm2 

21. 

Scattered  Pulse  Rate,  gen. 

V  =  n  A  v 

8  s 

250  sec'1 

22. 

"  "  »  det. 

vd  =  vg  “v/vmax 

5  sec  1 

23. 

Power  at  Focus 

PF  ’  WtWw 

Ps  .  Pp(d/2a  )z 

p  =  p  /16  r2 
c  s 

2.5  x  10"2  w 

5  x  10  !t  w 

24. 

Peak  Scattered  Power 

25. 

Collected  Scattered  Power 

2. 5  x  10  2  w 

26. 

Scattered  Power  to  Photo¬ 
multiplier 

PPM  "  'cWqVbWf 

-R 

1,4  x  10  w 

27. 

Scattered  Energy  to  Photo¬ 
multiplier 

Es  =  PPM  TS 

,  ,  *16 
4x1'  J 

28. 

Number  Scattered  Photons 
to  Photomultiplier 

Ns  =  Eg/hv 

1000 

29. 

Number  Scattered  Photoelec¬ 
trons  Generated 

Np  =  T]  Ns 

120 
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(19)  This  is  the  maximum  pulse  length  for  a  particle  passing  through 
the  axis  of  the  focal  volume., 

(22)  Here  v  is  the  maximum  velocity  corresponding  to  one  iree 
max 

spectral  range.  For  _\v  =  2GHZ  ,  vmay  =  570  m/sec  . 

(24,25)  The  scattering  cross-section  is  taken  as  the  projected  area  of 

the  (spherical)  particle  and  the  scattering  is  assumed  isotropic. 

The  aspect  in  which  the  actual  performance  is  most  likely  to  fall 

short  of  that  calculated  is  probably  the  focal  spot  size.  This  is  due 

imperfections  in  the  optical  components.  If  the  linear  dimensions  of  the 

focal  volume  are  in  practice  a  factor  a  larger  than  calculated,  then  the 

various  performance  fact<  rs  scale  in  the  following  ways.  The  spatial 

resolution  Ky  ,  the  projected  length  ,  and  the  scattered  pulse  length 

-g  increase  as  a  .  The  projected  area  Ag  ,  and  the  pulse  rates  vd 
2 

increase  as  a  >  while  the  scattered  powers  Pg  ,  Pc  ,  Pp^  decrease  as 
-2 

y  •  The  scattered  energy  Eg  and  the  numbers  of  photons  Ng  and  photo¬ 
electrons  Np  decrease  as  o” 

4.2.5  Alignment  and  Preliminary  Testing 

(a)  Initial  Testing  of  Laser  and  Fabry-Pernt 

As  supplied,  the  laser  operated  multi- line  and  multi- mode  and  it  was 
necessary  to  fit  the  prism  and  etalon  supplied  to  convert  it  to  single¬ 
mode  operation.  The  prism  is  inserted  into  the  optical  cavity  so  that 
operation  of  the  mirror  vertical  tilt  control  allows  one  to  select  one  of 
the  five  lines  in  the  blue-green.  The  etalon  is  also  inserted  into  the 
optical  cavity  and,  by  adjusting  its  vertical  and  horizontal  tilt  controls, 
a  single  axial  mode  may  be  selected.  The  mode  structure  is  illustrated  in 
Fig.  15. 
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Output 


Doppler  Broadened 
Line  Envelope 


Figure  15 

Laser  Mode  Structure 

The  exact  position  of  the  mode  pattern  within  the  line  envelope  depends  on 
the  length  of  the  cavity  and  thus  depends  on  the  mirror  settings  and  thermal 
tuning.  Four  adjustments,  two  for  mirror  tilt  and  tv-o  for  etalon  tilt 
must  be  made  to  obtain  maximum  output  on  the  strongest  mode  (nearest  to 
line  envelope  center).  Considerable  thermal  instability,  i.e.,  power 
variation,  extinction  anti  mode- hopping,  occurs  on  first  switching  on, 
but  after  a  prolonged  warm-up  (-.  2  hours)  the  laser  operates  stably  on  a 
single  mode  for  periods  of  1C  minutes  or  more  without  adjustment  of  the 
controls . 

The  Fabry-Perot  was  set  up  in  the  scanning  mode  with  a  period  T  = 

10  m  sec  and  synchronized  to  an  oscilloscope.  Initially  it  was  driven  from 
a  separate  ramp  generator  but,  because  of  problems  experienced  with  this, 
a  drive  from  the  oscilloscope  time-base  was  substituted,  using  the  circuit 
of  Fig. 16  to  control  the  amplitude  and  provide  a  d.c.  bias.  Use  of  the 
same  waveform  for  the  time-base  and  scan  drive  minimizes  effects  due  to 
nonlinear  waveforms. 

To  align  the  Fabry-Perot  and  test  the  frequency  stability  of  the  laser 
relative  to  it,  a  small  fraction  of  the  output  beam  was  directed  into  the 


>L  r«_ 
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Figure  16 

Fabry-Perot  Drive  Circuit 

Fabry- Perot  using  an  uncoated  optical  fie'  as  mirror  and  the  photodiode 
provided  with  the  Fabry- Perot  as  detector.  The  spacing  of  the  mirrors 
(2  GHz  etalon)  was  critically  adjusted,  and  the  orientation  of  the  inter¬ 
ferometer  relative  to  the  beam  set  to  give  the  sharpest  response.  A  half¬ 
width  of  <  10  MHz  corresponding  to  a  finesse  of  >  200  was  obtained. 

By  observing  the  steadiness  of  the  oscilloscope  display  with  the  scan 
amplitude  set  to  display  two  successive  Fabry-Perot  modes  (one  free  spec¬ 
tral  range)  as  in  Fig.  9,  it  was  possible  to  determine  the  frequency  stabil¬ 
ity  of  the  laser  relative  to  the  Fabry-Perot.  The  display  drifts  slowly 
due  to  thermal  turning  (primarily  of  the  laser)  and  exhibits  a  short-term 
jitter  at  a  few  H2  due  to  vibration  of  the  laser.  Once  the  system  is 
thoroughly  warmed  up,  the  thermal  drift  is  within  acceptable  bounds,  there 
being  negligible  drift  (i.e.,  less  than  the  resolution)  over  times  of  tens 
of  seconds  corresponding  to  the  maximum  storage  time  anticipated  (with  the 
slowest  pulse  acquisition  rate). 


The  short-term  jitter  causes  an  apparent  broadening  of  the  reference 
pulses  when  the  display  is  stored  for  times  —  1  sec,  compared  with  their 
width  with  single  scans  (10  m  sec) .  With  the  two  rc  Cerence  pulses  conve¬ 
niently  separated  by  8  cm  in  the  display,  the  width  increases  from  -0.5  mm 
on  single  scan  to  some  2  inn  on  storing  for  times  —  1  sec.  This  causes  a 
degradation  in  velocity  resolution  from  some  0.5%  to  some  2%.  While  this  is 
acceptable,  it  would  be  preferable  to  eliminate  it,  and  some  attempts  were 
made  to  do  so.  It  appears  that  the  laser  mirror  system  has  a  natural 
resonance  of  a  few  Hz,  which  is  excited  by  vibration.  Even  when  mounted 
on  a  very  stable  table,  so  that  externally  induced  vibrations  are  negligible, 
there  is  frequency  jitter  of  the  amount  noted  (-  40  MHz  total  excursion). 

This  is  apparently  excited  by  the  turbulence  in  the  cooling  water  channels. 

It  was  particularly  noted  that  if  the  water  flow  was  reduced  below  2.5 
gallons /min,  the  jitter  increased,  due  to  the  audible  onset  of  local  boiling. 
Apart  from  maintaining  a  flow  (3  gallons/min)  large  enough  to  prevent  this, 
it  seems  there  is  little  that  can  be  done  to  reduce  the  frequency  jitter, 
without  a  major  redesign  of  the  laser,  or  the  addition  of  an  elaborate 
feed-back  stabilizing  system.  The  mirror  resonance  at  a  few  Hz  is  also 
readily  excited  by  vibration  of  the  laser  mount,  as  could  be  observed  by 
lightly  tapping  the  stable  table.  In  comparison  the  Fabry- Perot  is  much 
less  sensitive  to  vibration. 

(b)  System  Alignment 

The  piocedure  devised  to  align  the  system  in  the  configuration  of 
Fig.  14  consists  of  the  following  basic  steps: 

(i)  To  make  the  beam  axis  coincident  with  u’e  axis  of  the  optical  bench, 

use  is  made  of  the  fact  that  the  beam  diame  er  almost  exactly  fills 


53 


the  telescope  entrance  pupil.  Thus,  by  sliding  the  telescope  carriage 
along  the  bench  and  watching  for  vignetting  of  the  transmitted  beam, 
one  has  a  sensitive  criterion  for  the  required  alignment  of  the  axes.: 
One  thus  adjusts  the  position  and  height  of  the  laser,  together  with 
the  position  and  tilt  micrometers  of  the  beam  splitter  to  obtain  no 
vignetting  as  the  telescope  is  slid  along  the  length  of  the  bench. 

(ii)  The  telescope  focus  adjustment  is  now  set  to  give  a  collimated  beam. 
This  may  be  done  with  sufficient  accuracy  at  the  present  stage  by 
making  the  beam  diameter  intercepted  on  a  distant  target  (~  5  meters) 
equal  to  the  diameter  at  the  exit  to  the  telescope. 

(iii)  The  next  step  is  to  obtain  a  reflected  beam  which  accurately  defines 

the  optic  axis  on  the  detector  side  of  the  beam  splitter,  allowing 

for  beam  displacement  by  refraction  at  the  beam  splitter  and  polarizer 

This  is  done  by  temporarily  setting  up  a  mirror  (an  uncoated  optical 
flat)  in  a  holder  beyond  the  telescope.  The  (visible)  beam  reflected 
back  through  the  telescope  is  split  by  the  beam  splitter,  and  the  mir¬ 
ror  mount  is  critically  adjusted  to  send  one  beam  exactly  back  on  the 

incident  laser  beam.  This  can  be  judged  very  sensitively  by  blowing 
smoke  into  the  beams  and  making  the  incident  and  returned  beams  coin¬ 
cident.  The  other  beam,  transmitted  by  the  beam  splitter  then  accu¬ 
rately  defines  the  axis  for  detection. 

(iv)  With  the  polarizer  in  position,  the  lateral  and  vertical  position  of 
the  Fabry-Perot  is  now  set  so  that  the  (visible)  beam  enters  the  en¬ 
trance  pupil  of  the  Fabry-Perot  centrally.  Using  the  photomultiplier 
as  detector  and  with  the  Fabry-Perot  scanning,  the  orientation  of  the 
latter  Is  adjusted  to  give  the  maximum  sharp  detected  signal.: 
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(v)  The  next  step  is  to  adjust  the  polarization  optics.  The  polarizer  is 
oriented  to  minimize  the  detected  signal,  under  Which  condition  it 
gives  maximum  rejection  for  the  vertically  polarized  laser  output,  and 

transmission  for  a  horizontally  polarized  signal.  The  quarter- 
wave  plate  is  then  mounted  to  the  front  of  the  telescope  and  oriented 
to  give  a  maximum  in  t  e  detected  signal.  The  reflected  beam,  having 
passed  twice  through  the  quarter  wave  plate  is  now  horizontally  po¬ 
larized. 

(vi)  The  spatial  filter  is  now  fitted  into  the  telescope  at  its  cross-over. 
Three  critical  adjustments,  two  for  lateral  position  ar  '  one  for  axial 
position,  are  made  to  put  the  pinhole  at  the  cross-over  of  the  inci¬ 
dent  beam,  as  Judged  by  the  maximum  uniform  illumination  in  the  trans¬ 
mitted  (expanded)  beam. 

(vii)  At  this  stage  the  cross-over  of  the  reflected  beam  will  not  lie  exactly 
at  the  pinhole  because  the  telescope  is  not  precisely  collimated.  This 
is  done  by  adjusting  the  telescope  focus  until  the  detected  signal  i. 
maximized . 

(viii)  The  temporary  mirror  is  now  removed.  The  signal  detected  at  this  stage 
is  small  and  constitutes  the  reference  signal  due  to  light  incidentally 
reflected  from  the  surfaces  of  the  beam  expanding  lens,  together  with 
light  scattered  from  the  edges  of  the  pinhole.  To  reduce  the  latter 
as  far  as  possible,  final  critical  adjustments  are  made  to  the  pinhole 
position  controls. 

(ix)  The  final  steps  involve  putting  the  focusing  lens  and  scan  mirror  in 
position  and  adjusting  them  to  direct  the  convergent  beam  centrally 
through  the  window  and  port.  This  can  be  judged  by  observing  through 
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the  MUD  channel  vith  the  exhaust  section  removed,  and  blowing  smoke 
into  the  duct.  The  telescope  focus  and  scan  mirror  vertical  tilt  are 
then  adjusted  to  bring  the  focal  spot  exactly  at  the  center  of  the 
electrode. 

When  the  scan  mirror  tilt  is  set  tc  put  the  focal  spot  exactly  on 
the  electrode  surface,  an  increase  in  the  detected  signal  is  observed. 
This  is  very  useful  since  it  gives  a  zero  reference  reading  on  the 
tilt  micrometer,  which  can  be  found  during  a  hot  run,  when  the  elec¬ 
trode  surface  may  be  expected  to  move  due  to  thermal  expansion  effects. 
The  magnitude  of  this  signal  scattered  from  the  electrode  depends 
strongly  on  the  state  of  the  surface.  With  a  smooth  shiny  surface  the 
signal  is  small  and  difficult  to  detect  as  an  increase  in  the  reference 
signal  due  to  incidental  reflection  and  backscatter  from  the  incident 
beam.  However,  with  a  rough,  oxidized  surface,  as  is  certainly  the 
case  in  a  hot  run,  the  reference  signal  increases  by  a  factor  of  2~-3 
and  gives  a  good  datum  for  determining  the  zero  position.  The  detection 
of  the  zero  position  from  the  increase  in  reference  signal  is  very 
sensitive  to  the  setting  of  the  scan  mirror  tilt  micrometer,  as  would 
be  expected  in  view  of  the  small  diameter  of  the  focal  volume.  The 
signal  peaks  in  a  range  of  micrometer  settings  about  -0.0001"  wide. 

Since  the  lever  arm  ratio  is  a  factor  —  7,  this  corresponds  to  a  ver¬ 
tical  deflection  of  the  focal  spot  of  less  than  0.001"  or  25  u.  While 
this  confirms  that  the  required  spatial  resolution  is  being  achieved, 
it  also  illustrates  the  sensitivity  of  the  boundary  layer  measurements 
to  vibration  of  the  mirror  mount  and  of  the  electrode  surface  during 
a  hot  run. 
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Calibration  of  the  vertical  deflection  of  the  focal  spot  in  terms 
of  the  mirror  tilt  micrometer  was  effected  by  placing  a  scale  vertical¬ 
ly  across  the  channel  and  observing  the  deflection  of  the  focal  spot 
for  various  micrometer  settings.  The  lever  arm  ratio  was  thus  deter¬ 
mined  to  be  a  factor  6.9. 

(c)  Preliminary  Cold  Tests 

A  number  of  performance  tests  were  made  with  the  system  set  up  on  a 

stable  table  in  a  dark  room. 

(i)  Focal  Volume  Dimensions:  The  diameter  of  the  focus  was  determined 

directly  by  measuring  the  waist  of  the  beam  from  the  side  using  a  high 
power  microscope.  By  blowing  smoke  into  the  beam,  the  diameter  of  the 
waist  could  be  observed  from  the  scattered  light  to  lie  in  the  range 
15-20  u.  This  is  a  factor  2-3  greater  than  the  calculated  value,  but 
is  entirely  adequate  for  our  purpose. 

The  measurement  of  the  intensity  distribution  along  the  length 
of  the  focal  volume  poses  a  difficult  problem.  Ideally  one  would  like 
to  suspend  a  single  micron- size  scatterer  in  the  focus  and  measure  the 
scattered  signal  as  the  particle  is  moved  along  the  axis.  Instead,  we 
mounted  a  rough  scattering  surface  normal  to  the  beam  on  a  translation 
stage,  and  attempted  to  measure  the  response  as  the  surface  was  moved 
along  the  axis.  However,  because  of  granularity  'n  the  surface,  the 
scattered  signal  was  quite  irreproducible .  Consequently  we  substitu¬ 
ted  an  optical  flat  for  the  rough  surface.  While  this  gave  reproduc¬ 
ible  signals,  it  is  clear  that  the  dependence  of  the  response  to  axial 
position  is  not  the  same  for  a  specularly  reflecting  target  as  for  a 
scattering  one.  Nevertheless,  such  measurements  give  some  indications 
of  the  effective  length  of  the  focal  volume. 
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Measurements  were  made  in  this  manner  with  a  number  of  different 


sized  pinholes  in  the  spatial  filter  to  determine  the  influence  of 
this  parameter,  as  shown  in  Fig.  17.  Theoretically  we  wild  expect 
the  axial  response  to  be  sharpest  for  the  smallest  pinhole  (4.7  u) 
but  in  practice,  better  results  were  obtained  with  the  largest  (15  |i) . 
Since  the  magnitude  of  the  reference  signal  scattered  from  the  edges 
of  the  pinhole  could  be  reduced  appreciably  by  using  the  largest  pin¬ 
hole,  this  was  adopted  as  standard.  The  length  of  the  focal  volume 
to  half  power,  measured  in  this  way,  is  ~0.8  mm,  which  is  some  five 
times  the  calculated  value  to  half  intensity  points. 

(ii)  Test  with  an  Incandescent  Electrode:  To  check  that  the  velocimeter 
would  allow  measurements  to  be  taken  adjacent  to  an  incandescent 
electrode,  the  beam  was  focused  at  grazing  incidence  on  the  surface 
of  a  tungsten  ribbon  lamp.  The  zero  position  could  be  found  by 

observing  the  increase  in  the  reference  signal  as  described  earlier. 

No  change  in  the  display  was  observed  when  the  lamp  was  heated  to 
~'2500°C,  which  shows  that  there  is  adequate  filtering  against  the 
background  radiation  even  when  the  spot  is  focused  on  the  electrode 
surface .; 

(iii)  Velocity  Resolution;  The  overall  velocity  resolution  of  the  system 
can  be  tested  by  placing  a  rotating  disc  at  the  focus.  However,  to 
obtain  useful  measurements  a  rather  high  disc  speed  —  100  m/sec  would 
be  required  and  such  a  test  was  not  made  since  a  suitable  motor-driven 
disc  was  not  readily  available.; 
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Relative  Response 


Influence  of  Pinhole  Size  on  Response 


(iv)  Tests  on  a  Water  Atomizer:  For  these  tests,  use  was  made  of  a  paint 


spray  atomizer  fed  with  water  and  powered  by  a  compressed  air  line 

at  50  psi.  The  nozzle  had  a  central  conical  needle  which  could  be 

adjusted  to  give  an  extremely  fine  spray,  such  that  the  water  consump- 

3 

tion  was  less  than  1  cm  /min  and  the  droplets  evaporated  within  a 
distance  of  some  10  cm  from  the  nozzle. 

The  beam  was  directed  at  an  angle  Q  —  26 0  to  the  jet  axis  and 
focused  to  a  point  about  1  mm  in  front  of  the  nozzle,  whereupon  strong 
doppler  shifted  scattered  signals  were  observed  in  the  oscilloscope 
display  as  shown  in  Figs . 18.  Signals  were  observed  over  a  wide  range 
of  velocities  from  —  50  m/sec  up  to  a  maximum  of  —  250  m/sec,  the  max¬ 
imum  being  set  by  the  velocity  of  sound,  since  the  nozzle  is  a  choked 
flow.  It  is  clear  that  the  flow  is  very  inhomogeneous,  there  being  a 
wide  range  of  particle  size  and  velocity  present,  as  well  as  a  high 
degree  of  spatial  inhomogeneity..  Although  the  characteristics  of  the 

flow,  i.e.  particle  concentration  and  size  distribution,  are  not 
sufficiently  known  to  allow  a  quantitative  evaluation  of  tie  veloci- 

meter,  the  atomizer  formed  a  convenient  source  to  prove  the  operation 
of  the  system,  at  least  qualitatively. 

The  conditions  for  the  results  of  Figs  18  are  as  follows:  Scan 
rate  1  msec /cm,  vmax  =  570  m/sec,  photomultiplier  voltage  2000, 
load  resistor  10  k  P  ,  vertical  sensitivity  50  mV/ cm.  Fig.  j.g  (a)  is 
for  a  single  scan  (10  msec),  while  Figs..  18  (b)  through  (e)  consist 
of  —  100  scans  corresponding  to  a  storage  time  —  1  sec.  The  narrow¬ 
ness  of  the  reference  and  scattered  signals  in  Fig.- 18  (a)  shows  that 

a  velocity  resolution  of  -  17»  is  achieved,  and  also  shows  that 

5  -  3 

the  particle  concentration  is  high  enough  (**  10  cm  )  to  give  of  the 
order  of  5  signals  per  scan,  The  degradation  in  velocity  resolution 
due  to  frequency  jitter  in  the  laser  is  evidenced  by  the  broader 
reference  lines  in  Figs.  18(b)  to  (e) .  The  only  difference  in  condi- 
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Fig.  18  Results  from  Cold 
Test  with  Water  Atomiser 


tions  in  these  figures,  apart  from  the  statistical  variation  from 
shot  to  shot,  is  the  setting  of  the  intensity  of  the  oscilloscope. 

This  illustrates  an  interesting  characteristic  of  this  type  of  display. 
Since  there  is  a  wide  range  of  particle  sizes  present,  with  more  small 
particles  than  large  ones,  there  is  a  corresponding  range  in  amplitude 
of  signals.  However  the  pulses  are  all  of  the  same  length,  determined 
by  the  time  constant  of  the  detector,  so  that  the  vertical  writing 
speed  for  large  signals  is  much  higher  than  for  small  signals.  Conse¬ 
quently  the  oscilloscope  preferentially  records  small  signals.  Fig. 
18(b)  shows  a  continuum  of  small  signals,  with  a  number  of  larger  sig¬ 
nals  showing  separately.  As  the  intensity  is  set  higher,  in  Figs. 

18(c)  through  (e) ,  si  the  infrequent  larger  signals  are  recorded. 

Figs.  18(d)  and  (e)  show  an  asymmetry  in  the  size- velocity  distribu- 

tion,  with  the  larger  particles  concentrated  at  lower  velocities. 

This  is  to  be  expected  since  large  particles  will  lag  further  behind 
the  gas  velocity  in  the  nozzle. 

It  is  clear  that  the  signal  to  noise  performance  in  Figs. 18  is 
entirely  adequate.  However,  in  the  absence  of  any  data  on  the  size 
distribution  of  water  droplets  present  in  the  atomizer  jet,  it  is  not 
possible  to  put  a  figure  on  the  signal  to  noise  performance  to  be  ex¬ 
pected  with  a  given  size  of  particle., 

(d)  First  Hot  Test  on  MHD  Flow 

As  a  first  trial  an  attempt  was  made  to  detect  signals  from  particles 
naturally  present  in  the  MHD  combustion  flow..  No  reliable  information 
was  available  concerning  such  particles,  but  since  the  MgO  firebrick  walls 
ablate  at  an  appreciable  rate  (~  1  mm/hour),  it  was  hoped  that  this  might 
provide  sufficient  scattering  centers.  In  the  event  no  signals  were  de- 
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tected,  and  subsequent  consideration  led  to  the  conclusion  that  since  the 
vapor  pressure  of  MgO  is  some  lOOtorrat  the  flow  temperature  (2400*0) , 
the  ablation  must  take  place  entirely  by  vaporization.  The  partial  pres¬ 
sure  of  MgO  near  the  channel  wall  probably  reaches  an  appreciable  fraction 
of  this  saturated  value,  otherwise  the  ablation  rate  would  be  much  higher 
than  observed.  The  possibility  of  condensation  of  MgO  to  form  particles 

near  the  cooler  electrode  is  unlikely  because  of  the  very  short  time  avail- 
-4 

able  (10  secs) . 

Thus  the  result  of  this  test  was  negative  and  we  are  forced  to  consider 
the  addition  of  particles  to  the  flow  as  discussed  below.  However  some 
useful  experience  and  information  was  obtained  in  other  respects. 

The  port  and  window  design  proved  adequate  and  no  window  damage  was 
detectable  after  approximately  one  hour  of  running.  The  laser  operated 
stably  with  regard  to  thermal  tuning  over  the  length  of  the  run.  It  was 
also  found  possible  to  follow  the  movement  of  the  electrode  surface 
(~0.005")  due  to  expansion,  by  monitoring  the  position  of  the  scan  mirror 
micrometer  to  give  an  increase  in  reference  signal,.  Some  indications  were 
also  obtained  concerning  the  severity  of  the  problems  associated  with 
vibration  induced  by  the  turbulent  flow.  These  are  twofold.  Firstly,  an 
appreciable  increase  in  frequency  jitter  in  the  display  was  observed,  and 
is  undoubtedly  due  to  vibration  communicated  to  the  laser,.-  This  degraded 
the  velocity  resolution  to  some  10%,  and  while  it  does  not  preclude  taking 
measurements,  it  is  clearly  undesirable.-  Probably  the  most  satisfactory 
solution  would  be  to  mount  the  laser  and  all  the  optics,  except  the  focus¬ 
ing  lens  and  scan  mirror  on  a  stable  platform  remote  from  the  MHD  rig. 
Secondly  it  was  cleai-  the  position  of  the  focal  spot  relative  to  the  elec¬ 
trode  surface  is  subject  to  appreciable  jitter,  which  would  degrade  the 
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spatial  resolution.  This  could  be  adduced  from  the  fact  that  the  increase 
in  reference  signal,  when  the  scan  mirror  micrometer  was  set  to  focus  the 
beam  at  the  electrode,  was  subject  to  jitter,  and  a  jittering  signal  could 
be  detected  over  a  wider  range  (-0.001")  of  the  micrometer,  corresponding 
to  a  jitter  amplitude  at  the  electrode  of  some  0.007".  This  may  be  at¬ 
tributed  to  vibration  of  the  electrode  surface  and,  or,  vibration  of  the 
mirror  mount.  The  main  contribution  is  probably  the  latter,  and  this 
indicates  the  need  for  a  more  rigid  mirror  mount. 

(e)  Seeding  the  Flow  with  Particles 

The  principal  requirements  to  be  met  by  the  seed  material  are  that 
it  should  be  readily  available  in  bulk  in  the  size  range  l-10u>  and 
that  its  vapor  pressure  at  the  flow  temperature  should  be  sufficiently  low 
that  the  particles  do  not  evaporate  during  their  time  of  transit  (—  1  msec) 
through  the  channel.  Of  the  various  refractory  powders  available,  classi¬ 
fied  for  size,  Alumina,  as  prepared  for  grinding  and  polishing,  appears  to 
be  the  most  suitable..  The  question  of  evaporation  may  be  assessed  as  fol¬ 
lows  . 

Assuming  that  the  evaporation  is  not  limited  by  heat  transfer  from 
the  gas,  so  that  the  particle  assumes  the  gas  temperature  T  (°K),  the 
mass  evaporation  rate  can  be  written 

E  -  2.7  x  1019(-^NJ  grams/cm2/sec  , 

where  p  (torr)  is  the  vapor  pressure,  m  (grams)  is  the  molecular  mass  and 
v.  =  1 . 38  x  10  ^  ergs/0  is  Boltzmann's  constant..  By  evaporation  the 
radius  of  a  spherical  particle  decreases  linearly  with  time  and  its  life¬ 
time  can  be  written 

E 
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secs  , 


where  p  (grams /cm  )  is  the  density  and  a  (cm)  is  the  initial  radius. 

For  Alumina,  we  have  the  molecular  weight  =  102  ,  p  =  3.7  and  at 
T  =  2700°K  ,  p  10  torr.  We  then  find  a  lifetime  t  =  34a  secs.  Thus 
for  lu  diameter  particles  the  lifetime  (1.7  msecs)  is  marginal,  but  for 
particles  of  5u  diameter  or  more  it  should  be  adequate. 

For  the  assumed  flow  conditions,  the  mass  flow  of  Al^O^  powder 

3  -  3 

necessary  to  give  a  particle  concentration  of  n  -  10  cm  ,  is  about  0.4 
grams/sec  or  4  lbs/hr.  There  would  appear  to  be  no  great  difficulty  in 
feeding  powder  into  either  the  alcohol  fuel  or  the  nitrogen  gas  at  this 
rate.  For  an  initial  experiment  it  appeared  to  be  rather  simpler  to  add  the 
powder  to  the  fuel,  and  a  mixing  system  was  devised  for  this  purpose, 
as  shown  in  Fig.  19  .  A  high  pressure  gas  cylinder  was  modified  to  contain 
the  charge  of  particles  and  fuel,  so  that  it  could  be  stirred  by  a  motor- 
driven  paddle.  The  mixture  is  driven  by  the  pressure  (150  psi)  of  a  nitro¬ 
gen  cylinder  through  a  flow  meter  and  mixes  with  the  alcohol  in  the  main 
fuel  line.  The  subsidiary  flow  of  particles  and  alcohol  constituted  a 
maximum  of  some  107„  of  the  total  fuel  flow,  so  that  the  particles  could  be 
added  without  changing  the  flow  condition  appreciably. 

(f)  Second  Hot  Test  (with  seeded  flow) 

For  this  test  the  mixing  cylinder  was  charged  with  10  gallons  of 
alcohol  and  10  lbs  of  Alumina  powder  of  a  nominal  5  micron  size..  At  the 


65 


To  Combuster 


Figure  19  Particle  Mixing  and  Injection  System 


maximum  flow  rate  of  --0.22  gal  /min  (-0.025  Ibs/sec  Alcohol)  indicated  by 

the  flow-meter,  the  particle  feed  system  would  run  for  approximately  an 

hour.  For  a  channel  flow  velocity  of  500  m/sec,  this  gives  a  calculated 

3  -3 

particle  concentration  of  —10  cm  ,  which  should  give  a  detected  signal 
acquisition  rate  of  a  f  per  second. 

Fig.  20{a)  shows  the  reference  signals  in  the  velocimeter  display 
before  starting  the  hot  flow.  Figs.  20  (b)  and  (c)  show  the  display 

after  warm  up,  with  the  flow  conditions  set  to  give  a  free  stream  velocity 

of  —  500  m/sec  at  a  cemperature  of  —  2400°C.  Fig.  20(b)  is  with  the  detec¬ 
tor  sensitivity  set  the  same  as  in  Fig.  20(a)  while  in  Fig.  20(c),  the 

amplified  gain  has  been  increased  by  a  factor  —  10  to  show  the  noise  level. 
These  figures  show,  from  the  width  of  the  reference  pulses,  the  degrada¬ 
tion  of  velocity  resolution  due  tc  the  frequency  jitter  induced  by  the 
vibration  coupled  from  the  M.H.D.  channel  to  the  laser. 

After  the  flow  had  become  steady,  the  particle  feed  was  acti  ited 
y 
y 

and  attempts^made  to  record  signals  due  to  particles.  However  no  positive 
results 'were  obtained,  and  it  was  thought  that  perhaps  the  particles  were 
being  evaporated  before  reaching  the  measurement  point.  Accordingly  the 
nitrogen  feed  to  the  channel  was  increased  to  reduce  the  flow  temperature 
somewhat  (to  approximately  2250°C),  but  still  no  signals  were  seen.  At 
this  stage  it  was  decided  to  try  adding  more  powder  to  the  mixing  cylinder 
which,  by  this  time  was  approximately  half  full.;  Since  no  more  Alumina 
was  at  hand,  5  lbs  of  fly  ash  graded  between  5  and  10  micron  was  added, 
but  still  no  unambiguious  signals  were  detected. 

At  this  point  it  was  thought  that,  since  the  nominal  flow  velocity 
was  500  m/sec,  while  the  free  spectral  range  of  the  Fabry-Perot  was  572  m/sec , 
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the  signals  were  being  «asked  behini  the  broadened  reference  signals. 
Accordingly  the  fuel.  Oxygen  and  Nitrogen  feed  rates  to  the  channel  were 
reduced  to  give  a  free-stream  velocity  of  some  250  ra/sec.  Immediately, 
doppler  shifted  signals  were  seen  as  shown  in  Figs.  20(d)  and  (e).  These 
were  taken  with  the  same  detector  sensitivity  as  for  Fig.  20(c)  and  cor¬ 
respond  to  storage  times  of  5  —  10  secs.  The  distance  of  the  focal  spot 
from  the  electrode  surface  was  10.4  and  6.0  mm  respectively  for  these  re¬ 
sults,  and  both  give  average  velocities  of  240  +  10  msec,  which  agrees 
well  with  the  calculated  free  stream  velocity  (250ra/sec). 

Some  comments  on  these  results  may  be  made  as  follows.  The  apparent 
velocity  spread  is  quite  large,  but  the  actual  velocity  spread  present  in 
the  flow  is  smaller,  being  artificially  broadened  as  a  result  of  the  fre¬ 
quency  jitter  due  to  vibration.  The  amplitude  of  the  signals  is  only  a 
few  times  the  noise  level.  It  should  be  noted  that  reducing  the  flow 
velocity  to  one  half  has  the  effect  of  increasing  the  scattered  pulse 
length  and  hence  the  detected  pulse  amplitude  by  a  factor  of  two.  The 
acquisition  rate  ofdetected  signals  was  also  lower  than  expected.  Both 
the  smallness  and  the  infrequency  of  the  signals  may  be  explained  by 
the  subsequent  observation  that  little  of  the  particulate  matter  was 
reaching  the  measurement  point. 

Following  these  measurements  in  the  free  stream,  measurements  were 
begun  in  the  boundary  layer.  However,  at  this  stage  the  run  was  termi¬ 
nated  due  to  the  appearance  of  large  and  destructive  lumps  of  solid  in  the 
flow.  Subsequent  examination  revealed  that  a  considerable  fraction  of  the 
particulate  matter  had  condensed  as  slag  on  the  cold  wall  of  the  combuster 
and  was  becoming  detached  and  projected  along  the  channel...  It  appears 
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chat  the  particles  were  melted  and  probably  vaporized  in  the  flame.  The 
centrifugal  action  of  the  oxygen  feed  then  threw  this  out  to  the  cold 
walls  where  it  condensed.  Unfortunately,  because  cf  the  sequence  of 
operations,  it  is  not  clear  whether  this  was  due  to  the  Alumina  or  the  fly 
ash  or  both.  However,  we  might  expect  that  the  fly  ash  contributed  largely 
to  this  since  its  composition  is  approximately  50%  SiC^.  25%  Al^O^  and  25% 
by  weight,  and  silica  boils  at  ~2200°C. 

Although  the  addition  of  particles  to  the  fuel  allowed  us  to  obtain 
some  preliminary  results,  it  is  clear  that  it  will  be  necessary  to  inject 
the  particles  into  the  nitrogen,  further  downstream,  where  they  will  not 
be  exposed  to  the  high  flame  temperature  or  the  cold  walls  of  the  corn- 
bus  ter... 

4.2.6  Summary  and  Conclusions 

Following  a  detailed  design  study,  a  laser  doppler  velocimeter  has 
been  constructed  for  measuring  the  gas dynamic  boundary  layer  velocity 
profile  adjacent  to  an  electrode  in  the  Stanford  combustion  MHD  facil¬ 
ity.  Preliminary  tests  indicate  that  the  design  requirements  with  respect 
to  spatial  resolution  (<  0.1  mm),  velocity  resolution  (<  2%  of  free- 
stream  value)  and  signal  detectability  (for  5u  particles)  can  be  met. 

In  the  first  hot  ^un,  no  scattered  signals  were  detected  and  it  was 
concluded  that  it  will  be  necessary  to  seed  the  flow  with  suitable  par¬ 
ticles.  In  a  second  hot  run,  seed  particles  were  added  to  the  alcohol 
fuel.  At  a  flow  velocity  of  500  m/sec  no  unambiguous  signals  were  ob¬ 
served.-  This  is  attributed  to  the  fact  that  they  were  only  marginally 
above  noise  and  masked  by  the  broadened  reference  signal  at  570  m/sec. 
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When  the  velocity  was  reduced  to  250  m/sec,  clear  signals,  several  tiroes 
the  noise  level  were  obtained.  The  velocity  deduced  from  these  agreed 
with  the  free  stream  velocity  calculated  from  the  known  flow  conditions  to 
within  the  experimental  error.  The  test  was  terminated  as  a  result  of 
problems  associated  with  the  formation  of  slag  from  the  particulate  mat¬ 
ter  in  the  combustor  section. 

While  these  results  are  encouraging  and  demonstrate  the  method  in 
principle,  it  is  clear  that  considerable  further  development  will  be  neces¬ 
sary  in  a  number  of  ar*as,  before  the  technique  can  be  reliably  used  to 
measure  velocity  profiles  with  the  required  accuracy.  The  aspects  in 
which  the  performance  must  be  improved  are  primarily  (i)  enhancement  of 
signal  to  noise  performance  and  reduction  of  the  reference  signal  amplitude, 
(ii)  improvement  in  the  frequency  stability  of  the  laser,  (iii)  reduction 
of  vibration  of  the  focal  spot  position  relative  to  the  electrode  surface. 
Furthermore  a  suitable  method  of  feeding  seed  particles  into  the  gas  flow, 
upstream  of  the  measurement  point,  must  be  developed..  There  is  every 
reason  to  believe  that  these  requirements  can  be  met  and  the  laser  doppler 
velocimeter  developed  into  a  precise  instrument  for  the  measurement  of 
boundary  layer  velocity  profiles  in  MHD  channels. 


5 . 0  Summary 


The  research  program  supported  under  contract  AF  F33615-69-C1171  has 
dealt  with  boundary  layer  studies,  with  current  distribution  and  leakage  studies 
and  with  boundary  layer  diagnostics.  In  the  first  area  a  computer  program  has 
been  written  which  couples  the  MHD  generator  core  flow  and  the  turbulent 
boundary  layer.  The  program  yields  the  local  heat  transfer  coefficient  and 
shear  stress  as  well  as  profiles  of  velocity  and  temperature.  Measurements  in 
an  argon  flow  showed  the  extreme  sensitivity  of  the  electrode  voltage  drop 
to  small  concentrations  of  CO^  and  N^.  An  extension  of  this  work  will  lead  to 
a  better  understanding  of  the  surface-sheath  component  of  the  voltage  drop 
which  has  not  yet  been  studied  in  detail.  The  influence  of  the  boundary  layer 
temperature  and  electrode  temperature  on  electrode  voltage  losses  for  several 
electrode  sizes  was  studied.  It  was  shown  that  under  some  conditions  small 
electrodes  had  less  voltage  drop  than  large  electrodes.  Designers  must  consider 
temperatures  and  spacing  of  both  insulators  and  electrodes  in  optimizing 
electrode  designs. 

The  current  pattern  in  a  conducting-wall  generator  showed  considerable 
current  collected  by  the  side  walls.  It  was  also  shown  that  Faraday  generator 
electrodes  which  extended  up  the  side  walls  resulted  in  more  power  than  conventional 
flat  electrodes.  A  calculation  of  current  patterns  showed  good  agreement  with 
experimental  data  at  the  anode  but  power  agreement  at  the  cathode  where  surface- 
sheath  voltage  drops  are  greater. 

A  new  method  which  uses  Faraday  splitting  of  laser  modes  is  suggested  for 
the  measurement  of  electron  concentration  in  the  boundary  layer.  A  laser  doppler 
velocimeter  and  its  application  to  boundary  layer  velocity  measurements  is 
described. 
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